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1¥ 


ry kK: 


an 
ml 
m1 
art 


design Mectpume Une Successful utilization of DDAM 
ta evaluate equipment response ta a design shock input ie 
dependent an the ability af the design shack spectrum to 
reflect accurately the structural environment existing in & 
given installation. Thies will be influenced Gy the size 
and weight of the submarine test subject, the Strie tur = | 
system emplayed in the coanmetructian af the submarine, and 
the inmteractian between the equipment, foundatian and hull 
structure used ta derive the design ehack spectrum. 
B. THE STATE OF CURRENT SHOCK DESIGN ANALYSIS 

Since its adaption, DDAM has been used withaut 
modificatian ta quality ship and submarine equipment 
inetallatiaon desiqnes. The original desiqn shack spectrum 
@eete 2) remains intact without revision. Shortly atter 
the adaptian at DDAM as & design requirement, vunderwater 
explosive testing af muclear submarines was Banned due ta 


the inherent riske invalved., Thies ban wae rel 


cr 


~, 
~ 


xed im } 7 


[Ts 
ay 


mn 


4 


for the low level explosive test af a N 462: 


Jase 


I) 


C 
Mm 


eybmarine. Mechanical and explasive teste af equipment 
inetallatians an shock simulation platforms are the anlx 
reqular check of desiqns meeting the standards imposed ty 
DDAM. A danger exists with this Kind of verification in 


that the tests may be desiqned ta accurately repraduce th 


Tr 


Input mations provided for in the design shock spectrum 


utilized with DDAM. 
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In recent years the increasing availability af large 
computers and the rapid development of numerical methads 
have pravided the engineer with mew tools ta use in the 
analysis af submarine shock response problems. 
Finite-Element/Finite-Difference methods allow for the 
analysis af structure and fluid responses and the work af 
Geers (Ret. 3] provides a means ta incarporate fluid- 
structure interaction effects. Several computer cades have 
been develaped ta analyze installed equipment response ta 


shock waves based an these farmulations. One at them is the 





Elastic Shock CELSHOK) cade developed by Weidlinger 
Assaciates under Detense Nuclear Agency and Office of Naval 
Research funding [CRef.4]). The mativation for this work has 
been to investigate modern submarine underwater explasive 
shack response im comjunction with a testing program using 
emall toa large scale models and shaped, explosive, tapered 
Charges. Throuqhcut its develcpment, ELSHOK has been 
Validated using these highly cantralled tests. In the low 
level explacive test af the SEN 456 class submarine, ELSHOK 
was used tao predict the level of equipment responses priar 
ta the tests and ta increase the level of caonmfidence 
associated with the tests. This cade has now been made 
available to the Naval Pastgqraduate Schaal by the Defence 


fluclear Agency along with support fram Weidlinger Assaciat 


tT 
if 


Peer RP OSE FUR THIS INVESTIGATION 


Peel] eee = Ghee |e l t= Sine ith ict provide a 


Ww 


qd 


Oy 
ate 
a 
a 


numerical method ber which the engineer can check fi 


“rr 
a 


ta ite adequacy for inatallation in a submarine. I 


i 
Th 


imtended ta be canservative and straightforward ta apply. 
When DDAM was first propased (Ref.4d], its authors warned 


that the desiaqn shack factors epeciftied (Ret.2) were nat 


ru 


explasive 


Ww 


absalutely determined and should be reviewed & 


ting progressed and data was accumulated. Although 


cr 
iT 
in 


te 


or 


if 


ting procedures for installed equipment components have 
become yverw sanphisticated, these simulators may mat accur- 
ately reflect the respanse of installations in modern da 


ewamearines. Without the abilitx ta carry aut high- 


imtenscity shack t 


7 


sta an real submarines, verification at 
the desiqn sheck epectrum used in DDAM can nat be 
accomplished as intended. The evalutian of numerical cades 
Capable af simulating the response of installed equipment 
ta a specified shock laading makes available a means ta 
check the applicability of the originally proposed design 
shack facters provided in DDAM to modern submarines. There 
Pene flexibility Built inte DDAM ta allow desiqn fer other 
than the originally prescribed shack input magnitude. CDOAM 


makes ma separate account for equipment“structure inter- 


The purgose of this investiqatian 12 ta examine the 


t 


response af several equipment madels inmetalled in & present 


dav seale, 47806 long ton (LT) submarine useing COAM and 
Saat ali The equipment madels range in weight from 


1,888-20 ,89800 pounds. Simple madele are used ta make 


mM 


arrelations possifble between the two methods and worst 


case results are abtained. Additionally, equipment“hul ) 


interactians are investigated for their poesible influence 


on averall equipment response. 
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Il. OYNAMIC DESIGN ANAL YTSIS METHOD 


A. DESCRIPTION OF DDAM 
As mentioned in the intraductian, DDAM is a simplitied 
modal analwsisa method which utilizes shack inputs which 


: oft 


ist 


were empirically derived from underwater explasian test 
realistic ship and submarine installations. Althaugh the 


amalysic if identical far both trpes of vesse] 


weakereter— 


ences im this document will] be to submarines. In order toa 


] assumptions must 


Ts 
je 


utilize DODAM to evaluate & desiqn, sever 


be satisfied, 


1.) The equipment and foundation make up A linearly 
elastic system. 


2.) The structure can be represented reasanabh] » 
well using = lumped parameter madel. 


3.) The structure is hard mounted ta the rest af 

the submarine. No batteming, or damped mounts 
are allowed. 

In Havy shock requirements, equipment shack classi- 
fications are based aon & Graded system. These grades run 
from A& for mission critical iteme, ta C far items that 
could indirectly affect the ability of the vessel] ta 


Functian as imtended. For instance, a lacker which becomes 


adrift and could injure persenne! ar surrpeliwenines aie ine a 


uf 


could be grade Cf. The grade of the item is epecified by 


contract and determines the deqree af shock qualification 

it must underga ta be canscidered sate far inetallatian. Tn 

thie study, all equipment models were considered grade A. 
The ideal way ta compute the elastic shack responce af 


an installed equipment is ta conzezider the entire structure 


etic eyrctem. The normal 


nu 


LSutmarine and equipment) as an e@] 


m can then be 


or 
a 


mades and natural frequencies af this 


= 


Getermined and the reepanse af each normal mode ta the 
water applied pressure laading= computed. The modal 
responses can then be superimposed ta qet the resultant 


eguipment response. In an everydar design methad such as 


f 


DDAM, this approach if impractical. Ta simplify the abi 
pracedure, ODAM utilizes the cancept at a shack epectrum. 
Anew point an the complete structure can be desiqnated as & 


reference point which can then be considered ta be a fixed 


Ck sist The denamic 


‘Th 


bace for the equipment an ane sid 
response af this portion can then be computed by 
superpasitian of normal made responses ta the mation otf the 
Fixed base. In DDAM, thie motion is specified as a =tep 
input tea each normal made and is calculated by empirical 
equations based on underwater explasian tests of madels and 
cubmarines (Ref. 3). A& detailed treatment af the 


mathematical basis for GCDAM is given br Butt CRef. SJ. 


In order to use DDAM, a lumped parameter model at the 
equipment muet be formulated ta atbtain the matrix 


expressions 


Ch) O03 + CR) Cee = 28 eae? 


we 
‘= l,' 
a] 


which hase the assaciated eiqenvalue problem 


ee a WS. = EK) 1X] a a? 


where CX.) is the matrix of eigenvectors and [w=] is the 


diagonal matrix af natural frequencies. The advanta 


ye 
Ww 
T 
a) 
“~* 


this procedure i¢ that it allows uncoupling of the 
equations af mation and subsequent ealutiaon of the denamic 
problem in terme of its component modes, Careful) 
formulation of the equations of matian ta capture all the 
Gominant characteristics af the equipment leads toa goad 
results with only a emall mumber af modes. When ODAM was 
developed, it was intended to be used primarily in hand 
eempeutations. The major difficulty encountered in the use 
of DDAM was the salutioan of the eiqenvalue problem which 
can readily Be handled taday with small computers. The 
nermal mode theory upon which DOAM is based js wel] 
established in shock and vibration practice. The area af 


cencernm if the adequacy of the adopted desiqn inputs toa 


in 


reflect present day submarine building practice 
Eee OeoiGbh SHOCK INPUTS 
The specification of design shock inputs for use in 


DDAM calculations determines the usefulnees of thie methad 


~~ 
~J 


in the evaluation af equipment installation designs. ODAM 
ie a standard. The desiqn input values must pravide 


caonmsietent fixed-based equipment excitatians tar the bBroaad 


A 
~t- 
it 
che 


base af submarine platforms in use today ta eimulat 


"standard" underwater explasgioaon. Ta simulate thi 


requifes 


ment here, charge weights and standoffs were selected ¢a that 


ifs 


each case corresponded ta a canstant value af energy Flux 
deneitywy 1.¢., 
WeRe = Conetant Ca 


In this equation, Wis the equivalent weight of the charge 


uC 


' 


in pounds af THT and F is the etandet? gisctance Fran eae 


ut 
ul 


hull in feet. The submarine platforms in use taodar ranqe In 


— “st 


size trem about 4,460-Ts.7 


i) 
i! 


ie Each submarine hull)] te 


J. 


characterized by its cqun modal properties mass dietri- 
bution, matural trequencies, mode shapes. The same piece 


mat equipment as is installed in & small submarine will 


demanstrate a different responce when inetalled in the same 


€Oti fl aij b x hai 


ou 


large suomarine. The effect af hull’ 
equipment interaction becames an important cansideratioan in 
large pieces af equipment tuned to ane or more matural fre- 
Quencies at the hull structure. The design shack input used 
with ODAM must incearporate these factors sa that when it 2 


said that a piece of equipment meets the specification for 


qualification, the qualification level is the same far all 


iT 


submarine classes, 


Leconte ies recearading WE, 1ittl= i¢ said 
met the Origins af the particular ehoack input epectrum 
Utilized ta evaluate equipment respanse., It hae een 
empirically derived ta provide values cansistent with the 


[mn the event that 1t 1 


if i 


ut 
iT 
ty 
wn 
m 
ae 


Seta Upon which it i+ 
Gecided ta design submarines ta resist underwater explasians 


mameocitterent sheek intensity than the one cheasen tor the 


~ 


present methad, a majar eftaort would Be required te 


canmetruct a data base upon which toa base the new ingut 
Walues, 


I. USING ODAM 


Ti 


Since the development af DDAM, new tools have rep] ace: 


the slide rules af engineers. Amang them are the readily 


fi 
y 
ce 
at 
es 


available desktop camputers, Cie tea e ee cet th lS Wick 84 t 


Was mecessary ta uee DDAM ta analyse simple structures art 


1Te 


fewer than 36 deqrees af freedam. Ta this end, & Computer 


i 


1 


v0 
ip 
't 


program has been written in the BASIC computer Janqua 
carry Gut the required computations. A listing and shart 


Users’ manual can be found in Appendix A. The pragram was 


ny 


veriftied by comparing results with hand calculatians and 


published sample probleme, 


TIi. EXPLOSIVE SHOCK RESP DNS ei i alae se Glee ccna 


&. GENERAL FRINCIFLES OF OFERATI ON 


“4 


he ELSHOK computer cade (CRetf. 4] calculates the 


trameient re 


if 


Te 


ponse af a submerged, ring-stiffened shell at 
revolution, with er without internal structure, tG an 


Underwater shack wave emanating from an explosive source 


or 


placed at an arbitrary location away from the ehell. The 


structure if cansidered ta be linearly laeti 


(Ts 
La 


=, and the 
surrounding fluid is treated as am imfimite acoustic 
medium. Medal structural analysis ifs used in all phases of 
the calculations. Interns&l] equipment response is treated 
by coupling the free-ftree modes af the empty ring-setiffened 
shell and the tixed-base mades of the interns] equipment 
through use of Gynamic Boundary canditians (Rets. /¢,d]. 

The etructure-fluid interaction is approximated by means of 
the Doubly Asymotatic Appraximation (CAA) due ta Geers 
CRef., 3]. The farm of the DAA used is that abtained when 
the normal fluid displacement mf the structure-fluid 
interface is expanded using surface expansion functians 
which are arthaqanal over the wet surface af the submarine. 
At frequencies of zero and infinity, the pressure-yvelacity 
relations are exact sa that in transient analysis, the DAA 


yields exact results at early and late times, and by the 


Zu 


narure aft it 


Wi 


Trimet henepreavl ges aremoath transt ti an 
between these two limite. Im effect, the DAA allawe for 
Pease ling of the fluid field pretlem trom the structural 
meelo protilem. 


t 


a! 
ay 


(Te 


ine =tructiralonpeblem salved By ELSHOK is cep 


au 
fs 


r: 
imto two parts. A modal substructuring procedure is used ta 
ealve the dynamic response problem for imternal equipment. 
The advantage of this is ta eliminate the need to handle the 
mades and matural frequencies af the combined estructural 
oroablem, at well ae the requirement for a combined srstem 


etiftfness matrix. Interactian f 


a 


rmees and moments, ano 
Compatibility of deformation at the shell-substructure 
attachment points are used ta salve far the dynamic response 
of the component parts. 

Reterring ta the tdeal case af DDAM tormulated wi thaut 
me use af a shock input value; ELSHOK is 4 numerical means 
to arrive at the input ta the Pine eacure wt thout depend- 
ence on an explasive testing database and with the added 
advantage that interaction effects between the hull and 


Sulestructure are taken imto account. ELSHUOK performs 


uy 


transient analysis whereas DDAM utilizes a simplified 
"front end" ta arrive at the maximum furces and deflectianse 
In & Qiven respanse problem for a single magnitude af 


loading. 


RB. ORGANIZATION AND IMPLEMENTATION OF ELSHOK 
The ELSHOK code is implemented a= %# feries af seven 


preagrams. The major components are: 


i BOSORS structural analwrzer far shell (Fet. 7] 


2. ACESNID - virtusl mast processor 
3. .PLPLASH = sivei’ + ie sae cee = Scr, 


structural analyrzer far substructure 
{Ref, j{&] 


ws | SPETCRUST substructure pracessoar 


é. USLOQOB - time inteqratian pracessar 


¢. , FPUSseasS 


platting processear 


Figure 1 depicts the general relatianship amang the 
code moadules. The numbers above correspond ta the 
chrancalaqical executian arder ta be fellowed in a given 


Malwsl 


0 
uf 
un 


In order to carry aut an analysis using ELSHOEK, 


a lumped parameter mode} af the submarine being analywzed 
must be formulated or available as well ae the equipment 
mode}. Meare will be said about madels later. The compu- 
tatians are carried out in four phases. 
B.i. Fhase I - Shel) and Fiuid Analysis 
BUSORG is a finite difference cade used ta 
capture the effects of the submarine hull) structure oan the 


overal}] response. Two models are actually used: a ful} 


model cantaining information about the major coanfiqurati oan 
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ha 
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ee ee as ae a as eee on am as ae ee 


(v) 
cy 


af the submarine and its general dynamic properties, and & 


Ts 


compartment model representing the Jlacalized ares oat 
Interest on the hull. The full madel is used ta capture 
the grass effects af the response such ase riqid Bady 
translation and whipping (HM = 1 madess, and railing and 
twisting “N = &@ torsianal madecs af the entire structure. 
The campartment model i¢ used to capture Jacal details at 
the response in the area of interest (NM = © Breathing and 
M22 mades>?. MM 15 the integral number at circumterential 


Waves ar harmonics in the circumferential dietributians aft 


the surtace expansian tunmctions. For each M, a 


wi 


eparat 


Tt 


BOCSOR4 calculation i¢ required. The BOSOR4G cade provides 
the in-vacuo free-tree modes and natural trequencies ot the 
submarine full. 


The secand part of phase To is the use at ACESHID 


in 
Cc 
Ww 
T 
CL 
= 
cr 
eats 
iT 


ta compute the virtual mass arrar, which i: 


lat time cantributiaon af the DAG. The virtual mase are 


‘Ti 
mt 
ag! 


ut 


determined from the solution, based oan 


Ww 


ImMple sources, 
of a steady-state low frequency problem in which normal 
displacements carresponding tao the surface expansion 


functions are applied ta the surface af & Cavity af 


‘Ts 


revolution in an infinite fluid having the same size and 
shape as the wet surface af the submarine. Unily one run aF 


this module is required for 41) values of MN considered. 


ie tinea ten ine nnace eet the “execution oat 
Beeson to compile the data from all BOSURG runs and 


ABACESHIO inte a shell-and-ftluid file tor later use. The 


w 


analxweis dane tar the wark in thi thesis considered anl» 


in 


the hull reétponece of ane submarine sa anls ane 
shell-and-fluid file was meeded for all subsequent 
calculatianese. 


B.2 Phase II - Subetructure Analwreis 





SAPIVU is a cammanly available finite element cade 


used in the ELSHOK cade ta perform the modal analyse 


mw 
Mm 
—- 


each substructure Cequipment moadel?d. The equipment 
inetallation is firet cast into an appropriate finite 
element representation and then SAPI'/ it used ta salve for 
the desired made shapes and frequencies. This intoarmatioan 
is placed in-a substructure made file along with the 
qeametry, stiffness and mats information. The masse, mode, 
and frequencyr informatian were also used in the DDAM 
analysic for ceamparative purposes. 


is the execution at the 


rye 
0) 
Ts 
— 
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The final step in ph 
PICRUST cade. This code prepares the input for phase IT! 
calculations by placing the above madal information on 4 
substructure file. Additianally, it calculates influence 
coefficients correspanding ta the base aor support matiane af 
the subetructure and accounts for the user-specified 


cannectivity. 


B.a. Fhase III - Time Inteaqratioan 





Following the completion at Phasee I and II, the 
chell-and-tiluid file is merged with the subetructure tile 
to form the input file +torpeeSskee which is Seon s) ea. 
inmteqration proececsar,. USLOB allows speciticatian at the 
charge weiqht and geometry. For the problems of cancern 


here, an exponentially decaying, empirical medel wae used 


ta represent the explosive shock source. I[t is af the 


+ Grit: 
PICR) = Kicw?7*/R9K2 exp (-t/9a) (4) 

where 3 

Pr¢R>? = the incident pressure at a radius R from 
the explasion can the hull? 

Ky =a multiplicative conetant fear incident 
pressure 

Kz = a spatial decay caonetant tor incident 
pressure 

t = time after arrival af shock wave at point 
Ot interest 

W =weiqght at epherical charge 

Bo = Ke Ww = 

Keg =multiplicative canstant tor time canstant 
of exponential decay 

Ka = 


spatial decar constant for time canstant 
oF exponential decay 


In all the cases, the geometry was such as ta maintain coan- 


vy 


etant energy tlux according to eq.¢3). The integration in 


time is done using a madified versian af the Runge-kKutts 
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mumerical imteqratioan echeme, During the integqratian, 


Bees CutpuUteSe Yeloacits histories for user-epecitied points 


em the hull and substructure. 


B.4. mba Sele eee ert t thi G 


Fhase IV capabilities were mot utilize in thie 


tigation. However, ELSHOK does have the capability ta 


meete Clot files econtaiming veleacity history information far 


prescessing on & Qraphics terminal. 


B.o. Detlectian Calculations 
As mentioned previously, ELSHUK produces velocity 


+ 


meer ies at the points af imterest om the substructure 


DDAM produces maximum deflections or forces at the same 


wi 
iyi 


paints for similar madele. In order ta make coampari¢soans 


Retween the two methods, it was decided to conmvert the 
meee ity histeries from ELSHUK tae deflectiam histories 


am .integrator based an Simpsoan’s 1/43 Rule. Far each 


equipment case, a Separate prag 
accamadate the variations im canfiguration between each 


madel and the ditterences im sutput fram the USLOR code 


ta point of interest specification. These programe are 


od 


,i¢ 


ediin Appendix EE. 
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ram was weitten im BASIL 


Wein 


ie tt 
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TV. MODELS USSR He vatieie eee 


ei. FINITE ELEMEMN FIN) PE Rtie gee ace ee 
The idea af reoresenting & given domain as a collectian 


ef discret 


LP 


elements ar values Wat mot first applied in 
fimite element/tinite difference methods. Early mathemati- 
Ciane estimated the value af T to nearly GH decimal places 
Ee pepresenting the circle with a polygon af a finitely 
large mumber af sides. The polrygan used by the mathemati - 
Cians if a madel of the circle ther wished ta represent. 
In modern mumerical methods, it 13 realized that very few 
pratlems in the real world can be analyzed without the wee 
of simplifying madele. BOSORG, the finite difterence cade 
im ELSHOR used ta describe the dynamic behavior at the sut- 
marcine hull, requires a greatly eimplified aca Ce le 
real world counterpart. Similariy, bath DDAM and ELSHOk 
Wee 4 substructure model] ta capture the fixed-base 
properties af installed equipment. 
RK. SUBMARINE HULL MODELING USING BOSOR4 

Actual submarine hulle are very complex etructures., 


However, ther do pottess certain characteristics that can 


or 
i 
+, 


be used ta simplite their modeling. From the aspec 
submarine hull shock response, the hull jf a free-free ring 


stiftened crlinder with suitable end clasures., The 


tJ 
CO 


Mmieerma!l Weights ef machiner>s and other cutfitting must be 
accounted far but these items can be smeared inta their 
adjacent hull structure. The madel thus becomes &« cylinder 
whase cross sections reflect the density af the submarine 
at a qiven frame. Since a cylinder is a surface af 
revolution, the properties of each cross section need only 
be specified at ane point. The reader will notice that any 
fee strijeture serving to stitfen the hull in a nen- 


axicymmetric sense, euch as nmon-circular frames aor massiy 


iT 


machinery foundations, can not be represented in thi 


un 


fashion. If these etructures are in the area of interest 
in the hull, they should be modeled as internal structures 
fame hase Il af ELSHOK. 

The BOSORG cade allows the submarine hull model to be 
constructed by specifying the properties of the submarine 
hull along «a line of revolution. A mumerical complication 
enters the problem when dealing with models of submarines 


because the structure aften repeats itself re 
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=A 
me 
NU 
a) 
= 
a 


its axis af revolution. With euch a canfiguration, ther 


Te 


re many modes in which the matian af the stiffeners jit oF 


ys 


small amplitude compared to that of the shell. The bavs 
between the rings vibrate at ftrequencies which may 
apporaximate thase carrespanding ta a simplyr-supported 
cylinder at the same geometry as the Bbayvw. Multiple or 


closelyw-spaced eigenvalues correspond tao modes in which 


Ni 
“NO 


one or more of the Bars 230 brat iio fiero ye eee 
unatftected. Toa avyaid eigenvalue difficulties, the orudent 


measure i¢ to analyse as small a segment at shell a 


he 


moesj ble in order ta avaid mumerical ditticulties with 
multiple, ar closely-spaced eigenvalues as indicated by 
Bushnell CRet. 7]. In submarine shack analysis, the 
imterest if in the grass response contributions at the 
entire hull ta a Given sectian af waupeiae ce The autheare at 
ELSHUK have allowed tor this through the use ot twa models. 
Be tull model is used ta abtain the grass contributians 
af the entire hull structure in the N = dttorsional? and 
f= | «translational and whipping?) mades. This will imclude 
such things as BulkKkheads, heavy ring stitteners, and end 
closures, A comeartment made | ie then used to abtain more 
detailed tntormation about the particular sectian at 
Imterest. Im this war, responses away fram the sectian af 
imterest do nat have to be carried through the remaining 


calculations since they do nat affect the response at this 


area anyway. 

Due to security restrictions placed in intormation fre- 
garding the construction techniques and arrangements ot U.5. 
muclear submarines, it was nat possible ta madel an actual 
submarine tor use in the analysis work contained in this 
thesis. With the aid at Weidlinger Assaciates oft New York, 
however, a madel resembling a general submarine hull struc- 


ture af S768 lang tons displacement was abtained. This 


3H 


model it available through the U.S. Naval FPostqraduate 
Penmeaci But 16 nat imeluded here sa that restrictions da not 
have ta be placed on the circulation at this wark. 

Sage wees so =Gnematic depletion ar the full subma- 
rine. The submarine being madeled has been cut inta 147 


slices aleng the tangitudinal axie. Im the model, intor- 


Mation regarding Young’s madulus, Paisson’s ratia, mas 


if 


density, thermal expansian coefficient, plating halt- 


thickness and whether or nat etitteners are eme 


re 


Teas ptra- 
wided far each ane of these points. Twenty-six discrete 
rings are represented along the axis of the model, twent> 
af which are in the local area of interest. These rings 
represent frames in the real submarine where T-stitteners 


Give increased stiffness to the hull. ifieSCeGtuality, these 


=i, 
ie 


etitteners are spaced irniy evenly along the entire madel, 
but smearing them outside the area otf interest will not 
aporeciably atftect the vibrational characteristics i thre 
area of interest. Using the BUSDORS cade, madal analysis 
using finite difterence techniques i¢ carried out. Table 

Ll is a listing af the natural trequencies and wave numbers 
tor each made retained in the analysis. Seven N = w 
torsianal modes, and thirty-eight N = 1 rigid bady 


dieplacement and whipping modes have been retained tar the 


fuj] madel. 
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ment tlecal area af interests model used in this anal» 
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In this model, all trames are madeled as di 
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Ts 


te Pine 


Stittener =. Since comoatibility must be maintained be tween 
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the twa madele, tcampartment and tull:, there | 4 One® te 


wi 
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Te 
if 


one carrespoandence in the area of interest between points 


im each model. The compartment model ist used ta obtain the 


madal cantributians corresponding ta N 2 2 normal, merid- 
lanal, and circumferential displacements). AS mentioned 


previqusl», deformations ar mations occuring beroand the 


a= 
ml 


region of interest should be unimpertant far the N 


O's" 


breathing modes and modes corresponding ta M 2. een 


bounding ring or bulkhead is very etiff in ite awn plane 
and maves primarily im that plane for lacalized laading. 


Table 2 is a listing of the matural frequencies and 


eircumterential wave numbers for ach mode retained in the 


Ts 


= 


analysizt. Far the compartment model, 25 mades carres- 
ponding ta N= 2, 24 modes carrespaonding ta N = 2, £4 mades 


carrespanding to N = 4 and 25 modest carrespoanding ta MN = ¢& 


vane 


have Been retained. A tatal af 254 modes were found Ey 
BUSOQR4 for the full and compartment modele. Of these 
mades, 1435 were retained for further calculations im phace 
Il and 111 of ELSHOK. The | Piemede= rif crear ances 


affect respanses primarily auteide af the area of interest. 
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Fiqure 4 i pees  Cieeaemen | ciel) Struc ture 


meee | and maotation canventians used in the remainder at the 
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ELSHURF analreiea. Two coordinate sretems eee, w,f ANG 


2 refer ta the submarine global coordinate trame arid “«.¥ 


8 


1 
sak 


Sie 2 Freter to the ubSstructure tactal cCeocrdinate trame. 


if 


both systems, the x-axes run from aft to forward and the 


Wt 


mean is Grientation is related to the 2 axis in the glabai 


in 


wetem Oy the angle “3 where @ refers toa the substructure. 


alt 
fe 


The variable = is used ta locate paints alang the meridian 
of the full medel, and u,v, and w are used for Jacal shell 
displacements and mations located circumtferentially By the 
angle 8. This cancludes the discussion af the submarine 
mode! used in ELSHUE. 

One af the conveniences realized in using =» ODAM 
th 


eke i absence of x 


analeyeis to carry out desi 


yw 
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mM 

Tr 
IT 


tT 


submarine model im the computations. The detigqn shack 


ea | U 


iTs 
in 


are intended ta provide all inputs to the substruc— 
ture model. The penalty which must Be weighed here is the 


moek er regard which thi 


int 


places am the peculiarities af 42 


iven submarine structure. On the ather hand, the ELSHUR 


WO 


suomarine model is difficult ta construct requiring a great 
dexl af information and skill te properly madel a qiven 


Rull, Gut the price is paid anly once for each submarine 


16 are 


rm 


lass. All subsequent calculations in this the: 


carried out utilizing the same shell model calculatians. 
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PeeeecUve ss RUCTURE MOCELING USING SAPIY 
Peter Vea sane onuh Calculatians, & substructure 


- 


model is formulated to determine the modest, mattez, and 


matural frequencies af the equipment eretem oeing design 


checked. When DDAM wat first developed, a computer program 


Was available which could salve the eigenvalue probiem for 


iT 


& 


cu 


up to twenty degr af freedom. However, the scarcity af 
camputer rescurces dictated that most uteri utilize hand 
computations far this purpose. Coansequentiy, mast earls 


CL 
‘Ts 
Ww 


mn checks were limited in scape. DDAM if well euited 


Wo 


to tabular methods af computation but with the common 
avVatllabillityw af small computers today, Rand caomputatians 
are no longer required, and the restrictian2 on degqrese at 
e | 


treedam are lar 


as 
~ 


removed. Thiet is mat ta say that great 


wo 


numbers af degrees af freedom are required in @#verw cate. 


uo 


Be caretul examination af the equipment ineatallation being 


checked, the majyor respanzse cantributions can be captured 
using & small number of modes. In the cases canmeidered in 


this thesis, emphasis wat placed an vweing equivalent model 


Let 


for DDAM and ELSHOR rather than closely modeling a real 
campanent in &@ submarine as it would be in each case in 


Cctual oractice. The interest here lies in fom the results 


he 


compare, far the same madels. 
Phase IT af ELSHOK utilizes the SAPIV finite element 


code to madel submarine-installed equipment. Because 


ELSHOK is written utilizing eioeqram megu te =e pene oe eee 


inmdependentl» and exchange data through cutput files, thre 


‘LI! 


SAPIU! code could be used ta run the medal analveize of 


equipment modele for ODAM and ELSHUOR thus ensuring equiual- 


ence im input tea both metheds. The models selected tar 
aAnNaRIKSIS were Kept emall ea that comparisons could be made 


ioe migra seive complication. 
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Three equipment models have been examined tar this 


thesis. Case IT] 16 a cantilever foundatian model with a 
1898 lbh weight attached ta its free end. Casez I] and Lil 


are a simple beam foundsatian suppearting twa weights at 
centerspan and having attachment paints on twee separate, 
Giscretering stitfeners. 

Cal. Subs truc tiir ese ase aed.) 


ected in 


ul 
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Ta 


The case | model develepment 1 
figure 3S. The mode! represents a 18908 lb valve supported Ey 
% cantilever foundation. The foundatian ic BOT t We yim 
AilSc CC 1x36 steel channels CRet. 11] and the valve is 
riqidly fixed ta the free end. Since the intent was not to 
Quality any particular value desiaqn, the entire valve has 


been modeled as 3 particle (lumped mzes >. The beam masts is 
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six lumped mastes. The foundation | 
desiqned to be fixed ta a discrete-ring stiffener in the 
submarine in an upright position. Although SJ-dimensianal 


motaian has been alleawed, the maximum deflection accure in 


the athwartship direction in response to a side-on shack 


AG 


Meeeing. ibmMiie “ase je Fepresentative of a reiativel»y Jigh 
equipment installation. 
if, Suyubetructure LCase Li Model 

The case [I model development if depicted in 
fiqure 4, The model reoresents a foundatian canetructed 
from AISC W 2Pxi?? I-beame (CRef. 11] supporting twa 18,88k 
lb weights. The foundation structure is represented fr 1s 
beam elements and 1? masses. Tt spans two discrete 
etiffenere to which ite ends are fixed. This mace] is 
representative af a large pump aor turbomachinery 
imetallation. 

fro, DouUbStmucture Case [Ii Made} 

The case TI] model] is similar toa the case TI] 
made] except the foundation structure has been changed ta 
AIlSt Lb] 2fxidd I-theams (Ref. 11] and the supperted weights 
have Been reduced to 3,888 Jb apiece. Thit made 
Pemresents an intermediate weight eguiagment inetallatian 
such as a main feed pump. 
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AISC C 186X386 Steel Foundation y 
Properties 


lee ascot ins! 
I = 206.0 in. 
YY 
J = 544.6 in.4 J 
web thickness 4.673 In. cr 


flange thickness = 6.9436 in. 
area of channel = 8.82 in.2 


Section A--A 


A A 
\ 
REAL WORLD FOUNDATION PINE SEEMEN aber 


AND MASS FOR DDAM & ELSHOK 


Figure 3 —- Case I Equipment and Finite Element Model 


AISC W 27X177? Steel Foundation CASE II 


Meee Properties J 
Tox 6748 oe 
yy 556 in.4 
J 20.1 in.4 x 
web thickness = @.725 in. 
flange thickness = 1.17@ in. 


_— 


area of beam S22 erin = 
Section A-A 
AISC W 27Ki1i14 Steel Foundation CASE ITI 


Properties y 
ley = 4898 inf 
Dis i59 in. 
YY 
J 7.36 in.4 x 
web thickness 6.578 in. 
flange thickness BZ S32 1s 


area of beam 33.6 in.2 
Section A-- 


FR 94 c=.) @ € 
——e FR 75 
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Figure 6 — Case II & III Finite Element Mode} 


Ve. ANSLYSIS 


The goal of the analysis procedure used im this ie 
tigation is two-tald. The primary chbjective We €G ama 
equivalent substructures using PPeiioand CLonBn -cl] aan 
parican af the maximum predicted deflectians can be made. 
secondly, any hull/“substructure interaction effects are ta 
be nated, The analysis 16 complicated by several tactare. 
All the calculations are based an constant energy flux which 
is dependent an charge weiqht and etandatf. Since DDAM 
shock tnputs ultimately represent the results ot explosive 


shack tests, the qeometry oat the “analytical charge" j 


un 


invisible ta the user and assumed in the empirical shack 
spectrum. ELSHOK uses inputs af charqe weight and standart 
ta calculate shock laading by eq. (4) sa the efttects aft 


their variation will change the transient response at the 


substructure even though the shock intensity is constant. 
In thie analysis, three sets af charge weight and standoff 
were used. é further complication results due to the fact 
that DDAM only can be used ta calculate maximum relative 
deflections or farces. ELSHOK calculates transient 
velocities incurred tyr the madel which must te canyverted ta 


maximum relative deftlectiane. 
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A. CASE ANALYSIS FROCEDURE 
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Pee eians be GUm Ee etifet cConetructing the 


equipment made | 


Wn 


Toe riives tigated Using seri. A dynamic 
analreia wad performed ta find the natural tréequenci#=: and 


made shapes of the madel. This information wat sav 


iT! 
CL 
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igs 


with the mass matrix for subeequent DDAM calculations. 
Using PICRUST, the SAPIY dats wat then reduced ta a suit- 
mee form ta be merged with the hull structure/tiuid data 
Meetmeoriate 1 calculations. Finally, all the resulting data 
Were integrated using the USLOB code ta compute the tran 
eremt Velocity pratiles tor each equipment installation 
comfiquration and charqe weightctandatff eet, 

The PICRUST code was used to specity the inetallation 
comtiquration of each model. Madel coordinate eretem orien-— 


tation toa the hull system } 


Ww 


Seer nietiny Cie aici e cc. 
fig. €49. The lacation of the equipment along the langi- 


Miai~mal axis is alsa 


pecified. Fear example, in the case | 


ml 


madel, the angle i 


Ww 


ef degrees and the base attachment i¢ 
at frame 93 in the submarine model within the compartment 
model). Far the case IT and ITI models, twa different 
arientatians were invesetiqated, ane with the model mounted 
athwartship and the other with the madel mounted verti- 
cally, i.¢., shock input from the side and battom af the 


submarine, 


a5 


The USLOB cade allows input at charge weight and 
etandott. Ta investtigate how these parameters attect the 
transient response at the substructure, three sets at val- 


ues were used tor these parameters based on three ditterent 


tit 


charqe weights. The parameters cha¢gen are listed in the 


following table. 


Charqe Weight S teidiewt s Dae tea nce 
(los TNT) flmches/Feet) 
a, 46 1,4149/117.7 
1H,Huu 2,400-"1466.7 
135,008 2,458/26049.2 


TABLE III - Charqe Weights and Standatts for Analyrsis 


Atter selection of all tactars attecting the gqeametry 
of the problem, the time step ine ment ane 1m Ceqg aan 


limite were specified. In each case, enaugh time steps 


in 


were chacen ta capture the peak response amplitudes. Thi 
number was tound by trial and errar. 


Atter calculation at the velocity pratiles for a given 


tal 


model, this intarmatian was integrated using Simpsan’s 1: 
rule ta obtain a deflection histary and the maximum detlec- 
tion response at the madel. As alluded to earlier, the 
models were canestructed ¢o that suitable ditterences taken 
between the velacity pratiles at designated points wauld 


* 


yield the relative deflection, at any instant in time, art 
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apeepaimt tm question. This precegure unfortunately re- 


Dea cePflenee tte et, Of sotimetry in the medel. For 


Mees, ta & 


Ww 


iietance, in the case I] model, the deflectian of the weight 
an the end af the cantilever in the athwartehip direction 
can De abtained By integrating the ditterence Between the 
velocity histories at the Base and the tip af the beam. 
Fiqures (Ff) and €8) show echematicall»y the differences 


taken ta calculate deflections for each model canfiqura- 


Vtip,athwartship 


SHOCK 
WAVE 


Vbase,athwartship 


Vrel = Vbase,athwartship - Vtip,athwartship 


Figure 7 — Relative Velocities Used to Calculate 
Relative Deflection Between Maes and 
Base for Case I 


V4,vert Vé,vert agent 
Va Ss 


SHOCK 
WAVE 


Vi,vert ¥ V3,vert V2,vert 


Vertical Relative Velocity of Nodes 3 & S wrt fixed base 
Nodes (shock wave from ATHWARTSHIP)D 


Va vert Vé,vert V5 ,vert 
V4 V5 
Vi,vert V3,vert V2 ,vert 
V4 Jo 


SHOCK WAVE 


Vertical Relative Velocity of Nodes 3 & 6 wrt fixed base 
Nodes (shock wave from below Keel) 


V3,rel = V3,vert — “i vert + V2, vert)/2 
Vé,rel = Vé,vert — (V4 vert + VS vert)/2 


Figure 8 - Case [IT & III Relative Velocity Calculations for 
Athwartship and Vertical Shock Waves 
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Perl was vsed ta calcul 
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e the maximum deflections using 
the infermation saved fram the SAPIV runs. A file wae con- 
Ssrpucted, using an editing program, to input the masses, 
matursa)] frequencies, and mode shapes for the madei ta the 
COAM program. During execution, the oragram utilizes 
prompts to guery the user about the twpe af installation 
and heck input directian. The maximum relative cGetftlectian 


iS calculated by the HRL tormula (Ref. i2]: 


CA = Sa aie x ~ 
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ii st modal deflectian of a 
an the mode} 


ctr=— 


‘elative detlectian at & point on 
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4, = maximum relative deflection of a point an 
the made] 


This formuls | 
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weed ta reflect the fact that ail mode=z da 
mot experience their maximum delections simultanecuel,. 
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configurations and shock input direction: 
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formatian used in the DODAM and ELSHOK calculations of thecte 


cases, Ae tull case |] analysis 3 
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presented in Appendix C. 
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I! Las I Inputs ta DOAM 


at Deqress at Freedom 
af Connected Degrees of Freedam 
ef Unconnmected Degrees at Freedam 


af Interface Paints 


of Frequencies In Analysis 


Substructure Natural Frequencies ¢cpes 


{ 445.41 al Guerre Se ee? oom 
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Although the mumber of cases and their permutations 


explored in this investigation are not numerous, ther have 


Seaem cantrived to accentuate the similarities and ditter- 
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ences between the GCOAM and ELSHORK methods af anal x 
mentioned previouely, ELSHOK provides a great deal more 


information sabeut the shock response af submarine-installed 


equipment than CDAM. However, ODAM it nat intended ta be & 
theeretical tool but rather is simply a design check. Far 
thie reasean, where ELSHDORK will predict different maximum 
deflections depending an charge weight and standatt, ODAM 
Will predict only ane value far all variatioane. 
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The predicted responses of DDAM and ELSHOK compare well 
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for this case. The ODAM calculation rielde 3 maximum de- 
#lectian af W.ALS3" for the end af the cantilever beam when 


subjected ta an athwartzehio shock input, whereas the ELSHORK 
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to Be the entire full mage, then the effective maze 


hull if sabeaut 418 slugs. The mass of the model is less 


than 2.792 | Wes ee eee eee 


ee ratia af abaut 8. Huh 


oy 


DOAM is based upen the appraximation that the huli mses can 
be considered infinite which implies a mate ratio af 
H.WUGBE... providing ma account far eeu ihe Gothen sone 
through the satete factor built inte its design spectrum. 
The clase agreement between ELSHOK and DDAM is probably due 


to the lack of response coupling in Oath methads af 


E. Cece 1] C2n, ane ste wre ac) 
Two equipment configurations were investigated in the 


ELSHUK analreesd far thi 
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case. In the problem qeametre in 
which the submarine i¢ side leaded, figure ¢7c2, the 
maximum deflectian is again in goed agreement with that 


predicted by DDAM; 6.2432 inches in the former versu: 


W.2794 in the latter. However, the input mation tea the 
structure is mot primarily in the directian of max imum 


deflection in this instance. In the secand configuration 
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Reere ction Wiss’. greater? is predicted By ELSHOK. There 
are two possible explanations. Same hullestructure 
interaction is toa be expected because the fundamental 
matural frequency af the substructure is 48.8 Hz which 
talle in the middle of the entire range of hull frequen 
eies. If, as in the former case, a hull effective mass of 
4168 slugs is ateumed, and the madel mass it 47.49 luge in 
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wpan which the shock input eoectrum used in DDAM | Cie: 


Most explosive teste are set up with the charge geometry at 
fiqure €Pc3 ta avoid problems AScUetatea ath Eieole migra— 
meen te the hull! due ta Baurancy erfects. Im the secand 
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factics Out perhaps mot simulated well br explosive test- 
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Since its adaption tr the maval shock community, GOeahH 


by 


ae provided an easy-to-use and convenient means ta desi gn 
check equipment inmetallations propeased for submarines. Tt 


iS independent af ane well defined submarine etructura! 


input and can be applied early in the averall design pra- 


cess Bbetore such details are well developed. However, the 
same Gualities that allew for the flexibilit Im DDAM ma: 


alsa cantribute toa possible imaccuracies when thie method 
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an all-encompassing evaluation af DDAM and ELSHOK. From 


the infinite mumber af possifle equipment configuration 


uw 


ate  Sh2 
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- three simple models have been selected which ex- 
amine aniy the full mounted equipment protlem. The intent 
af this investigation has teen ta examine the capability of 
the two methods ta predict the shack resnonse af these 


equipments and mat tao rank ane against the other.  ELSHQE 


sake 
rT 
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ar amy Varlatign et thi = cate =. muc complex and 


computer resource intensive to use ae x desiqn taal and was 


not intended tear this purpose, 


In the case af equipment whose response it relatively 
high frequency compared ta the hull response, DDAM and 
ELSHOK deflection predictians are in good agreement. The 
maximum hull response frequency is limited ta about Luu Hz 
Gene the &708 LT submarine used in this work. The gao0d coar- 
mebakion if attributed ta the lack oft interaction and 
coupling effects between the hull and substructure. When 
the response frequencies af the substructure approached 
those af the hull and the mass aft the substructure became 2 
siqnificant portion af the averall system, the coamparisans 
between ELSHUK and DDAM deflection predictions were mat 
qoad, 

The indications af this work are that DDAM does nat 
Serprectiy reflect the hull-substructure response coupling 


amplification which becomes apparent when the zubsetructure 
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response 1s tuned ta the hull respeanse. Additionally, the 
empirical database upon which DOAM is reliant for shack 


input may fail to represent the heavy equipment problem 


well, 


More basic research remains ta be done before many of 
the questions or criticisms posed br this thesis can be 
considered conclusive, Due to the scarcity of publiched 
ror ts describing the input data used ta generate the 
desiqn shack values used ty DDAM, the application af this 


method as a general design gualification shoyvld be reqgarde: 
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with caution. The Basie cp ie Desi an Sie ep = — een 
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CODAM should Ge available for analwrsis and revisioan. 
Additional wark remainge to be done ta determine where 


the discrepanci between ELSKOK and OGOAM tires 
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the substructure response soectrum Gut lend }i1 ttle | mraiee 
mation ta what happens in the transition region between the 


low-interactian and high-interaction frequency responce 
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However, eome madiftications are recammended here. The 
present Desiqn Shack inputs used im DDAM should be updated 


to reflect the ineresed 
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ize of modern day submarines. lf 


the details at the ariqinal formulation were KEncwn or 
revealed, the current values could Ge updated Gy generating 
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based an computer simulations. This 


new spectrum could then be yveritied by explasive testin 


Ww 


It if apparent from the hull-substructure coupling phe- 
nomenon nated here that a moditied shock spectrum shauld be 
formulated far each new class of submarine which would 
reflect the peculiarities of this clase versus the previcus 
Ones, 

At the present time, when new sutmarines cest billions 
of dollars and are by no means numerous, every effort must 


Fe made to ensure their survival in an underwater explasive 


shack enviranmennt. Thiet should tmelude taking & new Jaak 


at DOAM and ensuring it upholds the standards set by nmavai 
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The DDAM Program contained herein provides an ear 
means toa perform # small ta medium-scale DDAM analrvsis of a 


qiven internal equipment design. The program i¢ written in 


f 


TBP PO BASIC language and it nat restricted ta submarine 
anglwees amily Dut includes the standard DOAM case 
possibilities. The solution alqarithm faollowe the methad 
cited in Ret. 2. 
A. USING THE DDAM FROGRAM 
A.i. Optian Selectianse 
The DDAM program i¢ menu-driven and written ta be 
self-explanatory. Frier ta Starting he ei Geman. Fie 
be available an a diskette along with the BASICA. COM 
proar sm. 
Ta etart the program type: 
B 


AS 
FS 
F2 


4 
? 
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TLACENTER? 
*>DDAMCENTER > 
? 


The first command line laads and executes the BASIL inter— 
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preter and the secand line laads the DDAM pragram. The 
third line etarts execution at the DDAM Fragram. 
liponm starting, the title at the program ¢ODAM? 12 


wreeitten an the sereen and several secands later a meses 
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Pepe at = gesecribiing the previded tfunctians. Fallow the 
aicections. 
See ateeenmcp tiers Wilh appear far ceurce ct 


data and program exit to the operating system. They are: 


Ci] Input Mode Shapes and Mase Matrix tor the 
prablem fram Keyboard 


Peas entice shares ama tascs Platrix tram & 
Sst ile 


CS) iets less Hatrix ane akhitaness Matrix 
for previously tftormulated equatianse at matian 


bea Exit ta DOS 
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Mo optioen C4] if present; however, pravisians have been 
made ta inetall &# user-detined subroutine at line 22608 jn 
the program. 

Npetion CL) is usetul for amall problems and is mat 
recommended far esituations invelving mare than tour degrees 
of freedam. Mans al) datas ntl t (tGmibe carried ogy t 
using the keyboard. In larger problems, the entry af made 
shapes becames tediaus and error prone. 

Uptian C2) 16 Be ornie mee far most pratleme. fet 


allows the use of a disktile for input. The file i 
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weer-specitied and utilizes a free-format data structure. 
The intermation must tbe present in the file in a set order 


which is described later. 
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Selection of option Cs) causes aA message ta 
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Gispl as ‘tating that the JACOBI program must first te 

rum toa salve the eigenvalue problem for mode shapes and 
matural trequencies. This program lieting folloawe the ODAM 
listing im this appendix. IneéGrparating the JACOBI ar cee 
in DOAM would have reduced significantly the size af pron 
blem which could be handled. The JACOBI program is in- 


cluded at a convenience to the reader and wae not integral 


to the analrsis af the case studies done in this theel 


ui 


Tts operation if descrited later. 
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A.2: Uptionti) Exeeu ian 


Selecting aptian C1) allewe input of the macsee, 


ma 


modes, and matural frequencies af the prablem fram the 


kKerboard. The number at mades and masses are first input. 
A mass input subroutine queries the user interactively far 
the diagenal entries of the masse matrix. After input is 

completed, the diagonal elements are displared enabling the 


weer to make carrectiaonse as meceseary, mA zera carrectioan 
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terminates mass matrix input. Mext, the program asks for 


Input at the matural frequencies at the system. I 
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gGuires one nmatur. 


va 


1 frequency far each mode shape. The 
input format i¢ similar toa the one used ta input the masse 
elements. The matural frequencies are required in units of 
rPadians/ceecand. Each mode schape is then input. Atter 


initializing each mode, a COFPection Ghtian =sapaeiaea. 
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cer war cdtene at CxAlEWyates the participation factars and 
madal weiqhts far the protlem. 

Following these preliminary steps, OOAM requires 
Piesrmation ta fix the installation details for the equip- 


feat and trope at analysis to be carried gut. The standard 
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the 


Desiqn Wal 
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es are. calculated and ditplared as OAL)... DAC) 
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where m is the mumber at modes selected far the anal’: 


The Desiqn Values are the equivalent static acceleration 


calculated trom empirical formulas for each made. 


cr 


several autput options are provided. Madal tarces 


nf 


or deflections ar Death can be selected. After selectian, 
the Bpomorn | te calculatians are made and the results dis- 
pPlaved. Initially a print aption was separately pravided. 
However, im the present versian, selectian af the ditplar 
to screen aption with printer on will cause cutput tea qa tea 
the screen and the printing device. Disabling the printer 
Will cause output tea go only to the screen. Foallcawing 

the output options, choices are provided ta change the 
equipment inetallation txrpe and repeat the analywsizt aor tea 
exit the proqram. 
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Selecting oaptiaon [2] allows input af masses, mat- 
ural frequencies, and mode chapes from 2 disk file. After 
the data is read in, the program execution is similar toa 


that in option C1]. Option C2) prompts the user for the 


6? 


dimenztioan af the masse matrix, the mumber af madezs in the 


analysis, and the mame ot the input file. The input file 
is a free-format @SED) {1 le Sane ina tiene] = so > eon oe 


frequencies, and made ehapes for the problem. This tile 
may be caometructed from etructural analyzer CSePIV autput 
or any other problem formulation method using a suitable 
editer preqram., ‘CEDLIN, the editor provided with the eae 
Diek Operating System" is mat recommended tor thit purpace 
due ta its line arientation.) The format for the imput 


file follows: 


maeceCl> ...2. Maseem) (erie! 
emeqatio .... omeqatm) Krad’) 
Cede S leeetee) se cue aetna lenin 
Cet Gtiee bed? ea on Eee) 
Where m= number af masses and 


mn the analwsis 
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number oft mades 
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In the above format, the data is read in sequentially from 


left to right. The data it echoed to the screen .and“ar ., 
opeinter im order to provide a check for proper imput 


operatian., 
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B. PROGRAM LISTING - D DAA 


Westone y ORFs olate 6.1.8 
ZA PRIBT T FXXKERE REE FALH ESSE EKER E SS EEK K EERE RE ERA R ER ELLER EEA EELS RARER ED 
443° 
2) PRINT" +# 
+4" 
38 PRINT" #4 DDDD GOOD “As 
4%" 
48 PRINT" ** dD 0 0D 0 A A 
44" 
JO PRINT" £# dD 0 DD 0 AR 
43" 
68 PRINT" ¥# DD 60 DD =D A A 
£45 
78 PRINT" ¥# DODD DDDD A A 
xx! 
88 PRINT" xx 
*%" 
98 PRINT" +# DYNAMIC DESIGN ANALYSIS METHOD 
*3° 
188 PRINT" ¥% 
xx" 
118 PRINT" #4 
4%" 
{PB PRIRT 448% EH EEE EE EHR EK EEE ELE RRR EATER EERE SILER EEE EEE EASE EAE ERE 
444%" 
121 DIM MASS“ 48) 
122 DIM PHI<48,48) 
123 DIM WC48) 
124 DIM MX(48 48) ,Mx2¢48 48) SUM C48) ,SUN2°48) P4648) 
125 DIM WACd8) 
127 DIM YOC48) ,ABC48) , DA 48) 
128 DIM AAC4B) ,VAC4B) F448) ,DEFL 48) , DELTA 48,48) 
138 FOR 1 = | TO 5888 : NEXT | 
148 CLS: LOCATE 1,1,8 
{58 PRINT’ The purpose of this program is to allow the user a convenient meth 
od" :PRINT 
168 PRINT*by which to carry out a smal] to medium sized DOAM analysis of a give 
a" PRINT 
178 PRINT"piece of equipment. Familiarity with NORMAL MODE ANALYSIS and some u 
no" sPRINT 
188 PRINT"*derstanding of the theory behind DDAM is prerequisite to confident us 
e" PRINT 
198 PRINT"of this program. The program is prompt driven but the user should fay 
e "sPRINT 
288 PRINT"available to him / her either :°:PRINT:PRINT 


Mt 
rt 


a 


33 3 # 8 
= 


218 PRINT* f.) Mode Shapes and a Mass Matrix for the" 
228 PRINT* problem of interest.°:PRINT:PRINT 
238 PRINT® B.) AMass Matrix and Stiffness Matrix fron* 


7} 


248 PRINT? previously formulated EQUATIONS of MOTION’ :PRINT 
SPRINT: PRINT 


290 PRINT" ( Press ANY Key to continue )":LOCATE 25,1 ,8 
268 A$ = INKEY$; IF AS = ** THEN 268 
278 (ies 


288 LOCATE 3,1 ,8 

278 PRINT" In order to initialize the program input, choose one of the follow 
ing” sPRINT 

308 PRINT*menu selections Dd" SPRINT: PRINT 

318 LOCATE 18,28,8 :PRIN‘°{1] Input Mode Shapes and Mass Matrix’ 

328 LOCATE 11,28,8 :PRINT" for the problem from Keyboard’ 

325 LOCATE 13,28,8 :PRINT°{2] Input Mode Shapes and Mass Matrix" 

326 LOCATE 14,28,8 :PRINT® from a disk file ° 

338 LOCATE 16,28,8 :PRINT'(C3] Input a Mass Matrix and Stiffness Natrix’ 

348 LOCATE 17,28,8 :PRINT" for previously formulated equations of" 

338 LOCATE 18,20,8 :PRINT" motion® 

368 LOCATE 28,28,8 :sPRINT"(S) Exit to 005" 

378 A$ = INKEYS : IF AS d= CHRS649) AND AS <= CHRS(53) GOTC 488 

386 IF A$ = "* GOTO 378 

398 BEEP : GOTO 378 

406 8 = VALCAS) +: CLS 

418 ON 8 GOSUB 5888, 15888, 28088 , 25288 , 506 

308 SYSTEM :END 

SOQ 7 RRRKKAKRRE REESE EERE ERE SEER EE LES ESERIES IIIT 
681 ’ SUBROUTINE TO INPUT MASS MATRIX 

603 7333333 AGE IIS 
68) INPUT “Enter the dimension of the Mass Matrix ==) °; N 

618 IFN <= 28 GOTO 622 i 
628 PRINT:PRINT:PRINT"4*44 PROBLEM SIZE IS LIMITED TO 28 X 28 MATRICES +444 
* ":PRINT:PRINT: GOTO 688 

621 PRINT:PRINT:PRINT 

622 INPUT "Enter the number of modes you wish to use ==) "3M 

638° CLS: PRINT*#4##44 Enter the elements of the DIAGONAL Mass Matrix 4*4#4":PR] 
NT sPRINT | 

648 PRINT" If you make an error you will have an opportun-" 

638 PRINT" ity later to correct it.’ sPRINT:PRINT 

676 FOR 1 =1 TON 

686 PRINT "  MC*sEs","303") == “ssINPUT® ", MASSC(I) 

698 IF <] = 5) OR {7 = 18) GR {1 = 15) THEN CLS 

(86 IF MASS(i) <)> 8 GOTO 738 

718 PRINT’ 444% MASS MATRIX MUST BE POSITIVE DEFINITE --- INPUT NON-ZERO VALU 
E 444%" 

728 PRINT: GOTO 688 

738 PRINT: NEXT J 

748 CLS 

738 FOR 1 = i TON 

768 PRINT Aes oh iy oer eos 

778 NEXT | 

488 PRINT: PRINT :FRINT 

778 PRINT*Enter diagonal element number to correct mistaken entries or ° 

888 PRINT: INPUT " Enter < 2ERG > to continue ==)", RESPONSE 
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818 iF RESPONSE = @ THEN RETURN 

828 IF RESPONSE <= 8 GR RESPONSE > N THEN CLS : GOTO 758 ELSE | = RESPONSE 
898 PRINT:PRINT:INPUT “Enter correct element value ==) *, M¢J) 

846 PRINT:PRINT’MC*sT3?,":15") CHANGED TO °sMASSCT) 

8/70 FOR 1 = 1 T0 1508 : NEXT 1: CLS: GOTO 758 


888 RETURN 
[ROG “ RHKH KEKE HHH HHH KEKE EEE REE EEE EERE EERE EERE EEE EES EEE EEE EERE 
eee SUBROUTINE TQ INPUT HODE SHAPES 


[B28 ° RERHEE ARES RE REE ERECT EEE 
1958 CLS: LOCATE 5,1,8 

1866 PRINT" **##* Enter OMEGA’S corresponding to Mode Shapes #####£" PRINT 

1878 PRINT® f smaljest to largest )*:PRINT:PRINT 

1838 FOR 1} = 1 79M 

1698 PRINT ° omega¢*s13") == *ssINPUT® * Wd]) 

1188 IF ¢]=9) OR (]=18) OR (J=!5) THEN CLS 

1118 NEXT | 

PIZQeClsveruR t= 1 701 

1138 = PRINT"OMEGAC" 113") = Pslied), 

1148 NEXT ] 

1158 PRINT:PRINT sPRINT 

11468 PRINT "Enter frequency umber to correct mistaken entries or ":PRINT 
1178 INPUT ° Enter < ZERO > to continue ==>" ,RESPONSE 
1188 IF RESPONSE = @ THEN GOTO 1238 

1196 IF RESPONSE <= @ OR RESPONSE > N THEN CLS ;: GOTO 1128 ELSE I = RESPONSE 
1206 PRINT:PRINT:INPUT “Enter correct frequency value ==) ", WiI) 

1218 PRINT:PRINT*OQMEGA‘’:13") changed to ";W¢1) 

1228 FOR 1] = 1 70 2688 : NEXT 1 : GOTO 1128 

1238 CLS: LOCATE 3,1,8 

1248 PRINT**###* Enter EIGENVECTOR’s corresponding to Mode Shapes #*#**" sPRINT 
1298 PRINT" ( begin with Mode One )*:PRINT:PRINT 

1268 FOR 1 = 1 10H . . 

1276 =606FOR J =1 TON 


1288 PRINT" MODE SHAPE "323" Phic"yds") == "ssINPUT® * PHICI,J) 
1298 IF (J=5) OR (J=18) OR (J=15) THEN CLS : LOCATE 3,1,8 
1388000 NEXT J 


1385 CLS : PRINT:PRINT 

1318 PRINT "xesextteceteeee ODE NUMBER "sls" teneeeneeeeetee! 
1328 = PRINT: PRINT:FOR J = 1 TON 

1338 PRU ee = PRIN ee 

1346 = =NEXT J 

1358 = PRINT:PRINT 

1368 PRINT"Enter element number in mode to change or":PRINT 

1376 INPUT ® Enter < ZERO > to continue ===) ", RESPONSE 

1386 IF RESPONSE = @ THEN GOTO 1438 

1398 If RESPONSE <= 6 OR RESPONSE > N THEN CLS:GOTO 1318 ELSE J = RESPONSE 
1486 = PRINT: INPUT"Enter correct element value ==>"sPHI<] J) sPRINT 

1416 PRINT’Phit"sJ;") CHANGED TO *;PHI<I,J) 

1428 = FOR K = i TO 2008 : NEXT kK : CLS : GOTO 1318 

1438 CLS:NEXT | 

1448 RETURN 


MM TESISCSSSOSOSSSSOSSSISTICSCOSSSLIOCS ISS LOSS C COSC CS SSOSTOLOS COTS TTT T ELST TTT: 
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SAD “ «6 SUBROUTINE 78 DRIYE M0 bee es eee 

PUR TESIT LET I LILIES ITISI TTT TTITITELLLILeCETTSELTTTTSTEE Eee eet eeeeirirereteee: 
3628 LOCATE 3,1,8 

3838 PRINT" The program has been initialized to allow user input of the*:PRINT 
3848 PRINT"Mode Shapes and Mass Matrix. If you desire one of the other op-*:PR 
INT 

38598 PRINT*tions, PRESS < escape Key >; .....,0therwise ==)"sPRINT:PRINT 

3855 PRINT" PRESS < space bar > to continue )" 

3868 A$ = INKEYS : IF A$ = CHR3<27) GOTO 278 

3878 IF AS <> CHRS(32) GOTO 3868 

9875 CLS : LOCATE 5,1,8 : GOSUB 680 

3888 CLS : LOCATE 5,1,8 : GOSUB 1668 

3898 CLS : GOSUB 4888 

3188 CLS: LOCATE 25,1,8 

3118 PRINT 

3128 GOSUB 8868 

3138 CLS: GOSUB 8286 

3148 CLS : GOSUB 18888 

6888 7333333300333 
6818 ” SUBROUTINE 10 FIND PAR letra 

6828 ’ FACTORS 

6838 7 £33333 33a 
$848 LOCATE 25,1,6 

6858 7$(1) = TIMES 

6848 PRINT" CALCULATING PARTICIPATIGN FACTORS == start : “;TIHES 

6688 FOR 1 = 1 T04 

6698 = =8FOR J= 1 TON 


7088 MX(1,0) = MASSCJ) ¥ PHICI,J) 
7016 MX2¢1,0) = MXC1,d) * PHICI J) 
7628 = NEXT J 

7838 NEXT I 


7049 FOR I= 170M 

7058 = SUMI(1) = MXCI,1) 
7068 SUM2¢1) = MX2(I,1) 
7078 «= FOR J = 2 TON 


7088 SUMI4T) = SUMIAT) + MK(I WD 
1098 SUM2¢1) = SUM2(1> + MX2¢T dD 
7188 NEXT J 

7118 NEXT | 


7128 FORT = 170M 

7138 © PCI) = SUMIKT) / SUM2CT) 

7148 NEXT 1 

7158 FOR 1 = 1704 

7148  MBAR(I) = PCI) ¥ SUMICI) 

7178 NEXT I 

7188 LOCATE 25,48,8 : PRINT? done: ‘sTIMES 
7198 DUMMY = M 

7288 IF DUMMY ) 1@ THEN K/ = M/2 ELSE Ki = M 
7218 L=1:J8 

7215 LOCATE 5,1, 

7228 PRINT" #4244 PARTICIPATION FACTORS AND MODAL MASSES #4444" :PRINT :PRINT 
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Pez ecen Wiettt star Weber eeleteRo PND NODAL AAS SE S #448" 
7222 LPRINT 


7229 PRINT’ N = MODAL PARTICIPATION MODAL MASS * 
7226 PPS ="H# AH HHHR AHS” HH HHAR RHE OOS 
7227 LPRINT" N = MODAL PARTICIPATION MODAL MASS * 


7236 FOR } = L TO 

7248 © PRINT USING PPS + 1,PC1) MBARCI) 

7241 LPRINT USING PPS + 1,P(1) MBAR(I) 

7258 «(Jd t 

7268 NEXT | 

7270 IF DUMMY ) 18 AND J ¢ 11 THEN AS = INKEYS:GOTO 7298 ELSE GOTO 7311 
7288 = PRINT:PRINT:PRINT 

7281 LPRINT:LPRINT 

7298 = PRINT* PRESS ANY KEY TO CONTINUE...° 
7380 «IF A$ = ** THEN GOTO 7308 

7318 «KL=M: b= 18: CLS: GOTO 7215 


7311 PRINT:PRINT:PRINT* PRESS ANY KEY TO CONTINUE...° 
7328 A$ = INKEY$ : JF AS = "* THEN GOTO 7328 

7338 RETURN 

CABG % HEHE SEHK HH HFK EAS EAH ES BEEF EEE ES EAA EEE SESE EEE REESE EEE ESTES EASE 
8818 ” Soler nN Ueuetelnee | USCALUULATE NO bat 

8828 ” Weel b HTS 


BABA oH HHH EHH HHS S LEEK EAA ESAS SEES SSSR EERE EEE EERE EISELE EAST AE 
8835 CLS:LOCATE 3,1,8 

8848 PRINT"##4#4* CALCULATED MODAL WEIGHT S +exee" 

864) LPRINT"sx44%2 CALCULATED MODAL WEIGHT § #4444" 

8845 WAS = "HH HHH HHARHH” 

8846 PRINT : PRINT ° N WEIGHT IN XJPS* 

8847 LPRINT ° N WEIGHT IN KIPS* 

8858 PRINT 

8851 LPRINT 

8878 FOR 1 = 1 T04H 

8888 WACT) =386 * SUMiC]) # Pid> / 1088 

8898 PRINT USING WS 3 J ,WACT) 

889} LPRINT USING WS ; 1 ,WACT) 

8895 NEXT | 

8106 LPRINT:LPRINT:PRINT PRINT :PRINT" Press any Ke 
y to continue..." 

8118 A$ = INKEY$: IF A$ = "" THEN GOTO 8118 

8128 RETURN 

Q2AG (SHE HEEH HEE HAE HEHE E EEE E ESSE EEE EEE S EEE EERE EEE ESSERE EELS ESR 
8218 ” SUSROUL TNE tO ETE RHINE GEOMETRY 

8228 ’ OF ANALYS:iS AND AH ,V CH) , 0D TAD 

Q2IB HX EHEK KKH HAS SEEKER EELS SESE ESAS EERE EAE ESERIES SEER EEE 
8248 PRINT* Some information is now needed to complete the details about th 
o 

8298 PRINT’ type of vessel the structure J am analyzing is installed in and whet 
her” 

8248 PRINT*or not the analysis is to consider elastic or elastic-plastic defor 
ma-" 

8278 PRINT"tions." 
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8288 PRINT :PRINT 

8298 PRINT'PLEASE ENTER GNE OF THE FOLLOWING =>" 

8388 LOCATE 9,28,8: PRINT°{1) SUBMARINE \ hull mounted system )° 

8318 LOCATE 11,23,@:PRINT{2] SUBMARINE ( deck mounted system )° 

8328 LOCATE 13,28,8:PRINT"£3] SUBMARINE ( shel! plating mounted system )° 
8338 LOCATE 15,28,@:PRINT'(C4] SURFACE SHIP ¢ hull mounted system )° 

8348 LOCATE 17,28,8:PRINT°{5] SURFACE SHIP { deck mounted system )* 

8358 LOCATE 17,28, 9:PRINT'£C6] SURFACE SHIP ¢ she!] plating mounted system )°* 
8368 AS = INKEYS: IF A$>=CHR$(49) AND AS <= CHR$(54) GOTO 9398 

8378 IF A$ = "" THEN GOTO 9348 

8388 8EEP: GOTO 8368 

8398 B = VALCA$):CLS 

8488 PRINT*ENTER ===) [1] for ELASTIC ANALYSIS" :PRINT:PRINT 

8418 PRINT® [2] for ELASTIC - PLASTIC ANALYSIS ° 

8428 A$ = INKEYS : IF A$ d= CHRS(49) AND AS (= CHRS(58) GOTO 8458 

8438 IF A$ = °* THEN GOTO 8428 

8448 BEEP: GOTO 8426 

8458 C = VAL(AS) 

8468 IF C=2 AND (8 = 3 OR 8 = 6) THEN GOTO 2478 ELSE GOTO 8593 

8478 PRINT:PRINT:PRINT sPRINT® #***% ELASTIC - PLASTIC OPTION NOT AVAILABLE FOR #2 
¥*¥" sPRINT 

8488 PRINT" ###%% SHELL PLATING MOUNTED SYSTEMS HERES 

3431 LOCATE 18,26,8 

8482 PRINT" PRESS ANY KEY TO CONTINUE .... 

84383 AS = INKEYS : IF A$ = "* THEN GOTO 8483 

8498 CLS:60T0 8248 

8588 CLS:LOCATE 5,1,8 


8518 PRINT'ENTER =====) C1] for motion in VERTICAL direction’ :PRINT 
8528 PRINT {2} tor motion in ATHHARTSHIP direction® :PRINI 
8538 PRINT’ {3] for motion in FORE AND AFT direction’ 


8548 AS = INKEY$: IF A$ D= CHRS¢49) AND AS <= CHRS(51) GOTO 8578 
gju8 IF A$ = "" THEN GOTO 8548 
8568 BEEP: GOTO 3548 
35/8 0 = YVAL(AS) 
g571 IF B = 1 THEN LPRINT" SUBMARINE ( hull mounted system )"; 
8372 IF 8 = 2 THEN LPRINT" SUBMARINE ( deck mounted system )": . 
8573 IF 8 = 3 THEN LPRINT" SUBMARINE ‘ shell plating mounted system )°; 
$574 JF B = 4 THEN LPRINT" SURFACE SHIP ¢ hull mounted system )°; 
( 


8575 IF 8 = 5 THEN LPRINT® SURFACE SHIP ¢ deck mounted system )°: 

8576 IF B = § THEN LPRINT" SURFACE SHIP { shell plating mounted system )*; 
8577 IF € = 1 THEN LPRINT"® for ELASTIC ANALYSIS" 

8578 IF € = 2 THEN LPRINT® for ELASTIC - PLASTIC ANALYSIS" 

8579 IF 0 = 1 THEN LPRINTsLPRINT " for motion in VERTICAL direction’ :LPRINT 
8588 IF GO = 2 THEN LPRINT:LPRINT " for motion in ATHAARTSHIP direction" sLPRINT 
8581 IF 0 = 3 THEN LPRINT:LPRINT ° for motion in FORE AND AFT direction’ :LPRINT 


8578 FOR J = 1 TON 
8688 IF 8 <> 1 THEN GOTO 8688 


8618 ABCT) = 18.4% CC4BB+WACT) /(28+UATT))) 

8628 VOCT) = 28! & CCd8B+WACT))/ C1 BG+UACT))) 

8638 IF (D=1 OR D=2) AND C=1 THEN MULTI=1:MULT2=1 : GOTO 9118 ELSE 
8648 IF (D=3 AND C=1) THEN MULTI=.4:MULT2=.4 : GOTO 9118 ELSE 
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3658 
3648 
3478 
3688 
8698 
3788 
8718 
8728 
8738 
8748 
8758 
8768 
8778 
8788 
3798 
8888 
3818 
8328 
8838 
8848 
8858 
8868 
88/8 
8338 
3398 
8988 
8718 
8928 
8936 
8948 
8998 
3968 
8978 
3988 
3998 
9888 
7818 
9828 
9838 
9848 
9898 
7068 
9888 
9898 
9188 


IF <D=1 OR 02) AND C=Z THEN MULTI=1sMULT2=.5 + GOTO 9118 ELSE 
MULTIZ.4 + MULT2Z=.2 


GOTO 9118 
IF B {) 2 THEN GOTO 8788 


ABYT) = 3,2 * C4488 + WACTDI/C 28 + WACT))) 

VOUT) = 18! * CC4BB +aaCT) CIB B+WAT) )) 

JF (f=) AND C=1) THEN MULTI=1:MULT2Z=1 : GOTO 9118 ELSE 
IF (D=2 AND C=1) THEN MULTI=2:sMULT2=2 : GOTO 9118 ELSE 
IF (2=3 AND C=1) THEN MULTI=.8:MULT2=.6 +: GOTO 9118 ELSE 
IF (D=1 AND C=2) THEN MULTI=1sMULT2=.5 : GOTO 9118 ELSE 
IF (D=2 AND C=Z) THEN MULTI=2:MULT2=1! : GOTO 9118 ELSE 
MULTI = .8 : MULT2 = .4 


GOTO 7118 
IF 8 <> 3 THEN GOTO 8858 


ABUT) = 32 & CC4BRHWACT) /C284WACT))) 

VOC) = 188 * VC4BB HALT) /C1BB HALT) )) 

IF O=1 THEN MULTI=1:MULT2=1 : GOTO 9118 ELSE 
IF O=2 THEN MULTI=.2:MULT2=.2 : GOTO 9118 ELSE 
MULTI = .88 : MULT2 = .88 


GOTO 9118 
IF 8 <> 4 THEN GOTO 8758 


AGLI) = 28 * (437.5 + ACT) #C12 + ACT D/16 + WACTI*2)) 
VOC) = 68 * (412 + WACI))/(6 + WACT))) 

IF (D=1 AND C=1) THEN MULTI=1:MULT2=1 : GOTO 9118 ELSE 
IF (D=2 AND C=1) THEN MULTI=.4:MULT2=.4 + GOTO 9118 ELSE 
IF (D3 AND C=1) THEN MULTI=.2:MULT2=.2 : GOTO 9118 ELSE 
TF (DF) AND C=2) THEN MULTI=1:MULT2=.5 : GOTO 9118 ELSE 
IF <D=2 AND C=2) THEN MULTI=.4:MULT2=.2 +: GOTO 9118 ELSE 
MULT! = .2 : MULT2 = .1 


GOTO 9118 
IF 8 <> 3 THEN GOTO 9858 


ABCT) = 18 * 4037.5 + WACT)) C12 + WACT))/¢5 + WACI)*2)) 
Vacl) = 38 *# (C12 + WACTD/46 + WACI))) 

IF (D=1 AND C=1) THEN MULTI=1:MULT2=1 : GOTO 7118 ELSE 
IF <De2 AND C=1) THEN MULTI=.4:MULT2=.4 + GOTO 9118 ELSE 
IF (D=3 AND C=1) THEN MULTI=.4:MULT2=.4 : GOTO 9118 ELSE 
IF (=i AND C=2) THEN MULTI=1:MULT2=.5 : GOTO 9118 ELSE 
IF (DR2 AND C=2) THEN MULTI=.4:MULT2=.2 : GOTO 9118 ELSE 
MULTI = .4 3 MULT2 = .2 


GOTO 9118 


AGB(I) = 48 * 6637.9 + WACT))#012+ UALTDD/C6 + WACI))*2) 
VOCT) = 128 * CC12 # WACTD/66 + WALT) 

JF O=1 THEN MULTi=1:MULT2=1 : GOTO 9118 ELSE 

TF D&2 THEN MULTI=.2:MULT2=.2 +: GOTO 9118 ELSE 

MULTI = .f : MULT2 = 1 


9118 VOC1) = MULT2 * VOC1) 3 ABST) = MULT! * ABUT) : NEXT J 
9128 FOR 1 = 1 TOM 


9138 
9148 
9198 
9168 


AAY]) = ABCI) * 386 

VACT) = VOCI) *# WC) 

IF ABSKAACI)) > ABSCVACT)) THEN DAC]) = YVACT) ELSE DAC]) = AACI) 
IF ABSYDACI)) < 2316 THEN DAI) = 2316 


7? 


9178 NEXT ] 

W71 tier uUlalbegen.e 

9172 FOR ] = 1 TO Ms PRINT "DAC"S33") = “SDACH, 

9173 NEXT | 

9174 PRINT :PRINT :PRINT® Press any Key to continue ...° 
9175 AS = INKEYS:1F AS = ** THEN GOTO 9175 ELSE CLS : LOCATE 5,1,8 

9176 FOR ] = { TOM: LPRINT : LPRINT "DAC*31]3°) = "3DACT): HEXT To: LPRINT 
9188 RETURN 

{QB RQ 7 556363 323 3 EEE as 
18818 ’ SUBROUT IN Ft ee Sere) tae 

16828 ” FOR E.S  O7Re Gee eee eden 

{QQ3Q 7 36336336333 3G EE EAE 
10848 LOCATE 5,1,8 

19858 PRINT’ Almost done.... Please ENTER the number of the type of informat 
on* 

18668 PRINT:PRINT* you desire as output ===)" 


18878 LOCATE 18,1,8 :PRINT® [1] MODAL FORCES" 

19830 LOCATE 12,1,8 :PRINT® (2] MODAL DEFLECTIONS® 
19898 LOCATE 14,1,8 :PRINT [3] FORCES & DEFLECTIONS 
19895 LOCATE 16,1,8 :PRINT® [4] CHANGE ANALYSIS TYPE? 
1989 LOCATE 18,1,8 :PRINT? [5] EXIT To 00s" 


19100 A’ = INKEY$: IF A$ d= CHRS(49) AND AS (= CHRSC53) THEN GOTO 10149 
18118 IF AS = ** THEN GOTO 1a10¢ 

10138 BEEP : GOTO 10108 

18148 FLAGI = VAL(AS) + CLS 

19150 IF FLAG! = 2 THEN GOTO 19220 

19155 IF FLAGI = 4 THEN GOTO 16148 

10156 IF FLAGI = 5 THEN SYSTEM 

19148 REM +x*4%% CALCULATE MODAL FORCES xxxxxax4 

19178 FOR] = 1704 

19188 FOR J=1T0N 

18198 F(},J) = PCI) * MASS(J) * PHICI,J) * DACI) 
19288 NEXT J 

19218 NEXT J 

19220 IF FLAG! = 1 THEN GOTO 12000 

18238 REM x¥43%% CALCULATE MODAL DEFLECTIONS sxxxs4exs 

1a24a FOR I=1704 

19250 FOR J=1T0N 

19268 DELTACI J) = PHICI,J) * PCI) * DALI) 7 WKI)42 
19298 NEXT J 

19308 NEXT 1 

19318 FOR J= 1 TON 

19328 BIGGEST = @ 

19338 FORI=1T04 

1934@ IF BIGGEST < ABS(DELTA(I,J)) THEN BIGGEST = ABS(DELTA(] ,J)): Kel 
19358 NEXT J 

19348 S$ =a 

19378 FORT =1 70M 

19380 IF 1 = K THEN GOTO 104a0 

19399 = § + DELTACI,J)*2 - 

19408 NEXT J 


8 


16418 DEFLCJ) = ABSSDELTALK,J)) + SQR(S) 
18428 NexT J 

12888 REM *4444444% OLTPUT RESULTS *4444 444444 
12018 CLS : LOCATE 3,1,8 


12028 PRINT *° ENTER [1} TO ECHO RESULTS TO SCREEN®:P 
RINT : PRINT 

12836 PRINT * [2] TO ECHO RESULTS TO PRINTER®: 
PRINT :PRINT 


12648 AS = INKEY3 

{2858 IF A$ d= CHR$(49) AND AS <= CHR$(50) THEN GOTO 12078 

12668 IF A$ = "* THEN GOTO 12848 ELSE BEEP: GOTO 12848 

12878 FLAG2 = VAL(AS) 

12088 IF FLAG2 = 2 THEN GOTO 12308 

12898 IF FLAG! = 2 THEN GOTO 12198 

12188 FOR ] = 1 TOM 

iit Chom Oeale o,1,8 3 PRINT 449 NODAL FURTES FOR MODE 
Cr els*) Clbs) *4*" sPRINT 

WAMMGUPRINVe EER: #420 OAL FORCES FOR HODE « *sI;* ) [lb 
s] #42" : LPRINT 

12112 LPRINT " N ": LPRINT 

12114 PRINT ° N° : PRINT 

{2115 FF = "HH «HAL RRRHHHCC**® 


i2120 FOR J=1 TON 

12138 PRINT USING FFS ; J,FCI,J) 

12131 LPRINT USING FFS ; J,FCI,J) 

12148 NEXT J 

12158 PRINT :PRINT sPRINT :PRINT 

12168 PRINT * Press any Key to continue ...° 
12178 A$ = INKEYS: IF A$ = "* THEN GOTO 12178 

12188 NEXT | 


12198 IF FLAG] = 1 THEN GOTO 12289 

12288 CLS: LOCATE 5,1,8 

12218 PRINT’ #44 DEFLECTIONS AT MASS ATTACHHENT ?0Q 
TN 7 S #4#" PRINT: PRINT 

12211 LPRINT : LPRINT"## DEFLECTIONS AT MASS ATTACHHNE 
NT POINTS ##°:LPRINT 

12215 003 = "WH RHRHHH 

12228 FOR 1 = 1 TON 

12238 PRINT'DELTA ("315°") = °3 

12231 LPRINT" PEETaee 13") = "> 

12232 LPRINT USING DDS ; DEFL(I) ; 

12233 LPRINT*in." 

12235 PRINT USING DDS ; DEFL<I) ; 

12236 PRINT “in.* 

12248 NEXT | 

12258 PRINT :PRINT:PRINT 

12251 LPRINT:LPRINT 

12268 PRINT’ Press any Key to continue ..,' 
12278 AS = INKEY$ : IF A$ = "" THEN GOTO 12276 

12288 CLS : LOCATE 5,1,8 

12298 GOTO 16668 


(9 





12386 PRINT ° OPTION NOT AVAILABLE YET;===) USE <¢ Prtsc > KEY ‘ 

12318 FOR i = 1 TO 3888 +: NEXT 1 GOTO 12888 

LBBB % HERRERA HHS SHREK SEAR E SNEAK AEE ESE ERE RLS RRE LEER RRR EIS 
15610 SUBROUTINE 18 1 NPT sie 5S 7 fh dee 

15828 PROM D1°S kh oe ee 

L5Q3B o RH HHH HHH HHH HHH HEHEHE R EERE ELSES SESE E SAREE ES EES ESSERE AS EEE SS SSS 
15848 CLS : LOCATE 5,1,8 

15858 PRINT* The program has been initialized to allow user input of the’ sPRINT 
15868 PRINT°Mode Shapes,Frequencies, and Masses via an input file." :PRINT 

15878 PRINT : PRINT"I4 you desire one of the other options," :PRINT:PRINT 

15888 PRINT* PRESS (escape Key>: ......otherwise ==)":PRINT : PRINT 
15898 PRINT* PRESS <space bar> to continue ° 
15188 AS = INKEY3 +: IF A$ = CHR$(27) GOTO 278 

15118 IF AS <> CHRS¢32) GOTO 15188 

15128 CLS:LOCATE 5,1,8 

13138 INPUT “Enter the dimenston of the Mass Matrix ===> °3N 

13148 IF N <= 38 THEN GOTO 15178 

15158 PRINT:PRINT:PRINT**####% PROBLEM SIZE JS LIMITED 70 58 DEG OF FREEDOM *4%< 
3! 

15168 FOR ] = 1 TO 1888 : NEXT 1] : GOTO 15138 

15178 PRINT:PRINTSINPUT"Enter the number of modes you wish to use == 
15175 PRINT:PRINTSINPUT*Enter the input file name = 
191868 CLS 

15198 73033333333 EE IE 
19288 ’ INPUT MASS,FREG, AND PHI 

15218 73433333333 EEA 
15228 OPEN FINS FOR INPUT AS #1 

15248 FOR 1 = 1 TON 

15258 INPUT #1, MASS‘1) 

15268 NEXT | 

15278 FOR 1 = 1 7T04 

13288 = INPUT #1, Wod) 

19298 |= Wel) = WET) * 6.283185 

15388 NEXT | 

15318 FOR J=1 TON 

15338 «FOR ] = 11048 

13348 “INPUT #4 ,PHI<I,J) 

15358 NEXT | 

15348 NEXT J 

15378 733333333 
15388 ’ OUTPUT MASS, OMEGA, AND PHI 

15398 7 3333333333 EAB a a 
15488 LOCATE 5,1,8 : PRINT*MASS MATRIX DIAGONAL ELEMENTS * : PRINT 

13481 LPRINT"###44# HAS S DIAGONAL ELEMENT S #4x0" 

13482 LPRINT 

15485 MM$ = "HH HH. HHHHH**** . 

13418 FOR 1 = 1 7T0N 

15428 = PRINT USING M43 ; I ,MASS(I) 

13421  LPRINT USING MMS ; I ,MASSCI) 

15438 NEXT | 

15431 PRINT:PRINT* ooo PRESS ANY KEY TO CONTINUE® 


~ 


ya ie 
ye 
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15432 A$ = INKEY$ : IF AS = ** THEN GOTO 15432 

15434 CLS 

15443 PRINT:PRINT :PRIMT 

15458 PRINT*FREGUENCIES USED IN ANALYSIS RPS)" :PRINT 

15451 LPRINT 

15452 LPRINT'Hxeee FREQUENCIES USED IN ANALYSIS (RPS) 
444%" ° PRINT 

15460 FOR 1 = 1 TOM 

15478 PRINT 13" si) 

{5471 LPRINT Te? 9 "sW¢1) 


13488 NEXT | 

13481 PRINT:PRINT” oe ePRESS ANY KEY TO CONTINUE’ 
13482 AS = INKEY$ : IF AS = "* THEN GOTO 13482 

19484 CLS 


13498 PRINT: PRINT sPRINT 

13491 LPRINT:LPRINT 

19988 PRINT*MODES USED IN ANALYSIS" :PRINT 

{0081 LPRINT*##42# MODES USED IN ANAL YSIS ##¥8€"SLPRINT 
19018 K=M : US = "HEHHARH*SS** 

13528 IF K > 3 THEN GOTO 15688 

13538 FOR J=1 TON 


15548 FOR 1 = 10K 

15558 PRINT USING US; PHICI,J); 
1555! LPRINT USING US; PHI‘I,J): 
15568 NEXT 1 

15574 PRINT °° 

15571 LPRINT * * 

15588 NEXT J 


19598 = PRINTsPRINT 

13591 LPRINT : LPRINT 

13592 PRINT: PRINT" 1. PRESS ANY KEY TO CONTINUE’ 
13993 A$ = INKEYS : IF AS = "" THEN GOTO 15593 

13994 CLS: GOTO 15748 

$3688 =FOR Ji TON 


15618 FOR 1 =1705 

15628 PRINT USING US: PHICI,J): 

15621 LPRINT USING US: PHI<I,J): 

15438 NEXT | 

15648 PRINT * 

154d1 LPRINT * * 

15658 NEXT J 

15468 PRINT:PRINT:LPRINT :LPRINT 

15461 PRINT ve. ePRESS ANY KEY TO CONTINUE? 


13662 A$ = INKEY$ : IF AS = "* THEN GOTO 15662 
13663 CLS: PRINT: PRINT 
13678 =FOR J = 1 TON 


13688 FOR 1 = 6 TON 

19698 PRINT USING US; PHI<1,J); 
13691] LPRINT USING US$; PHI(1,J); 
19788 NEXT ] 

19718 PRINT * * 


8] 


SANs LPRINT ° * 


15728 = NEAT d 
15738 = PRINT :PRINT:LPRINT sLPRINT 
1373) PRINT® vee PRESS ANY KEY TO CONTINUE? 


15732 A3 = INKEY3 : IF A$ = *" THEN GOTO 15732 

15734 CLS: PRINT : PRINT 

15735 CLOSE #1! 

15748” RXR ELH EHA HERES ERASE EERE EEE EERE EASE RELEASE EEE ERE EE EERE EERE REE AE LEER 
15758 ° NORMALIZE MODE SHAPES 

15 75B 7 HEESRHS HRS ERAS ARES REESE RELA RE EEE ERE EERE EES EEEE EES EE EES SEER EERE EEE 
15778 FOR] = 1 104 


13788 = BIG = 8! 

19798 «=6FOR J=1 T0N 

13888 IF BIG ¢ PHIC],J) THEN BIG = PHI(],J) 
13818 NEXT J 

15828 «=6FOR J = 1 TON 

15838 PHI(] J) = PHICI,J)/B36 

19848 = =NEXT J 

13858 NEXT | 


LSB GB © RXR KHY SHEESH AE HH ESAS HAE ELL ES SEES EGGERS LEE EEE EE ER ES REESE ERE EES 
19878 ° OUTPUT NORMALIZED NODE SHAPES 

L5G SHAKES H AHS KEKE HERE KEKE EEK EERE EERE REE E REE E REESE EEE EEE EEE 
15898 PRINT'NORMALIZED MODES USED IN ANALYSIS" :PRINT sLPRINT 

15891 LPRINT"#ex4*4 HORMALIZED MODES FOR ANALY S 1] 5 4taxe" 
15892 LPRINT 

19988 K =¥ 

15918 IF K > 5 THEN GOTO 15988 

19928 FOR J=1 TON 


15938 FOR] = 170M 

15949 PRINT USING US; PHI<I,J); 

15941 LPRINT USING US:PHICI J); 

15958 NEXT I 

15968 PRINT * * 

1596! LPRINT © * 

15978 NEXT J 

15971 PRINT: PRINT" ve eePRESS ANY KEY TO CONTINUE" 
15972 A$ = INKEY$ : IF A$ = "" THEN GOTO 15972 

15973 CLS 


19988 «= LPRINT : LPRINT : GOTO 16128 
13998 »=6FOR J=1 TON 


16808 FOR 1 = 1705 

16818 PRINT USING US; PHI(I,J); 
16811 LPRINT USING U$; PHI<I,J); 
14828 NEXT | 

14838 PRINT * * 

16831 LPRINT * * 

16848 NEXT J 

16858 | PRINT:PRINT:LPRINT:LPRINT 

16851 PRINT? vs sePRESS ANY KEY TO CONTINUE" 
16852 A$ = INKEY$ : IF A$ = "* THEN GOTO 15852 
16853 CLS 
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16868 FOR J=1 T0N 


14878 FOR 1 = 4 TON 
14898 PRINT USING US: PHI(I,J); 

14081 LPRINT USING U3; PHICI,J); 

14898 NEXT | 

14188 PRINT * * 

1418! LPRINT * * 

16118 NEXT J 

16128  PRINT:PRINT:LPRINT:LPRINT 

16121 PRINT" ve PRESS ANY KEY TO CONTINUE! 


14122 A$ = INKEYS ; IF A$ = "* THEN GOTO 16122 

16138 CLS : GOSUB 4888 

16148 CLS : LOCATE 21,1,8 

14158 PRINT 

16155 CLS ; GOSUB 9888 

14168 CLS : GOSUB 8288 

14178 CLS : GOSUB 18808 

20888 CLS 

20018 LOCATE 5,1,8 

20028 PRINT'#4# YOU MUST FIRST USE THE PROGRAM JACOBI TO SOLVE EIGENVALUE +##" 
20038 PRINT"#4* PROBLEM. IT WILL CALCULATE NATURAL FREQUENCIES AND MODE #+#° 
20848 PRINT'#%* © SHAPES THAT CAN THEN BE USED TO RUN THE DOAM PROGRAM, #42" 
20858 PRINT: PRINT: PRINT: 


20868 PRINT® =n aes ENTER 4F3>JACOBI<RETURND® SPRINT 

28878 PRINT" {F2>" 

28888 STOP 

29088 CLS:PRINT**#% USER MAY USE OPTION 4 TO INSTALL HIS/HER SUBROUTINE ##3° 
23818 PRINT® 44% BEGINNING AT LINE 25088. ¢%%* 


29828 FOR 1 = 1 TO 3888 : NEXT 1: GOTO 278 
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C. PROGRAM JASOB] > THSt eae ales eine se 


The following BASIC program mar fe d 


~— 


j=. 


‘fl 
tT 
ws 
snl 
ft 


mT 
modified bye the user ta salve the Eiqenvalue problem ance 


the equatians of matian are faa la ter see 


im 
fi 
rte 


ee ee 
atory in execution with the &xcenpticgns that theo. mae 
replaces the stiffness matrix and the B matrix replaces the 
mass matrix. The program only requires that the Goeper Stn 
angie elements of the stiffness matrix be supplied alana 


with the disaqonszl elements af the mass matrix. 


D. PROGRAM JACOBI LISTING (Adapted from ref. 13) 


4 CLS 

D REM ==CYCLIC JACOB] METHOD= 
18 DIM AC18,18), 8418) ,UC18,18) 
19 REM ==INPUT DATA== 

26 GOSUB 78 

28 REM==INITIALIZE= 

36 GOSUB 668 

48 GOSUB 988 

49 REM ==QUTPUT RESULTS= 

06 GOSUB 758 

68 GOTO 3888 

65 REM ==INPUT DATAS= 

76 PRINT 

86 PRINT * CYCLIC JACOB] METHOD" 
168 PRINT 

188 PRINT 

196 PRINT " ENTER THE SIZE OF PROBLEM (i.e. no. of rows)" 
268 PRINT 


rR 
ti 


218 INPUT "N= "tN 

228 LPRINT “PROBLEM DIMENSICN = 's N 

238 INPUT "NUMBER OF SIGNIFICANT FIGURES = *; Si 
248 LPRINT 

250 PRINT 

268 PRINT * ENTER THE UPPER TRIANGULAR ELEMENTS OF MATRIX A ¢ by columns ?:° 
278 PRINT 

280 LPRINT "UPPER ELEMENTS OF MATRIX A =" 
298 LPRINT 

324 FOR J= 1 TON 

338 PRINT 

348 PRINT * ENTER UPPER PART OF COLUMN *; J 
358 PRINT 

$68 FORT = 1704 

378 PRINT © AChsTe* tege8) = fs 

388 INPUT ACI, J) 

381 LPRINT * ACTSTs* 88) = 8: ACTS) 

382 LPRINT 

398 NEXT I 

488 PRINT 

448 PRINT 

458 NEXT J 

460 FOR 1 = 1 TON 

478 = 1-1 

480 FOR J = 170 1 

498 ACI,J) = AGT) 

508 NEXT J 

S18 NEXT J 

520 PRINT 

530 PRINT * ENTER ELEMENTS OF DIAGONAL MATRIX B:! 
531 LPRINT " ELEMENTS OF DIAGONAL MATRIX B ° 
532 LPRINT 

548 PRINT 

578 FOR =170N 

580 PRINT * B('sTs*) = *: 

598 INPUT BCT) 

591 LPRINT * B¢'sIs") = "3BC(1) 

$08 NEXT I 

$18 PRINT 

611 LPRINT 

b58 RETURN 

$55 REM = INITIALIZE TOLERANCE AND MAX NO. OF ROTATIONS = 
68 2 = 2% SI 

478 T1 = 1/(18*2) 

688 PRINT 

$98 CLS: PRINT "TOLERANCE = 37! 

$91 LPRINT "TOLERANCE = ";Ti 

708 R= SN 

718 Ri = 8 

726 12 = 4! 

738 Ni=N-1 
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748 RETURN 

746 REM ==QUTPUT EIGENSOLUTIG = 

756 PRINT 

768 PRINT 

778 PRINT "SGLUTIGN’ 

771 LPRINT 

738 PRINT 

781 LPRINT 

782 LPRINT "x4 xxx ee xxeeeeE SOLUTION 40464445 eee Hee ER! 
783 LPRINT 

888 PRINT " NUMBER OF ROTATICNS REQUIRED = ":R1 
881 LPRINT “NUMBER OF ROTATIONS REQUIRED = “;R1 
882 LPRINT 

883 LPRINT 

818 PRINT 

828 PRINT 

838 FOR J = 1 TON 

848 PRINT 

858 PRINT " EIGENVALUE "sJs" 1S "sBcJ) 
891 LPRINT "EIGENVALUE “3J3" 1S * 384d) 
832 LPRINT 

853 LPRINT "ITS EIGENVECTOR IS ° 

868 PRINT 

878 PRINT " JTS EJGENVECTOR IS ° 

888 PRINT 

898 FOR J] = 1 TON 

988 PRINT UCI,J) 

981 LPRINT UCI J) 

982 LPRINT 

918 NEXT J 

928 PRINT 

921 LPRINT 

938 PRINT " PRESS SPACE SAR TO CONTINUE 
948 JF INKEYS = ° ° THEN 942 

941 GOTO 948 

942 CLS 

958 PRINT 

968 NEXT J 

978 RETURN 

973 REM ==E]GENPROBLEM SOLUTIONE= 

988 GOSUB 1138 

985 REM ==PERFORM ONE CYCLE OF ROTATIONS= 
998 GOSUB 1298 

995 REM == CHECK TOLERANCE = 

1888 JF X1 < T1 THEN 1118 

1985 REM == CHECK NO. OF ROTATIONS = 
1818 IF R1 > R THEN 1848 

1828 72 = .1 # X! 

1838 GOTO 998 

1848 PRINT 

1958 PRINT ° #£% ERROR ###° 
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1851 LPRINT "Set eese see eee eee eee ERROR 4HHKKHE REESE EEE EERE! 

1868 PRINT 

1878 PRINT ° NO COMVERGENCE ATTAINED’ 

1871 LPRINT “ NO CONVERGENCE ATTAINED ° 

1888 PRINT 

1881 LPRINT 

1898 PRINT * WITH ":Ri:* ROTATIGNS® 

1891 LPRINT * WITH °:Ris* ROTATIONS" 

1168 END 

1118 GOSUB 1838 

1128 RETURN 

1138 FOR 1 = 1 70N 

11448 FOR J=1 TON 

1136 UC1,J) = 8 

1158 NEXT J 

Me uC) eed 

1188 NEXT | 

1198 FOR 1 = 1 TON 

1288 B1 = S@R(BC1)) 

1218 B11) = 1 / 8! 

1228 NEXT | 

1238 FOR | 

1248 FOR J 
J 


1248 NEXT J 

1278 NEXT | 

1280 RETURN 

1298 Xi = @ 

1388 FOR K = 1 TON! 

1318 Kl = K+! 

132 FOR L = KI TON 

1338 Al = ACK,K) 

1348 A2 = ACK,L) 

1358 A3 = ACL,L) 

1348 X = A2 ¥ A2 /CAL # A3) 
1345 REM == CHECK IF ROTATION IS NEEDED = 
1378 IF X ) Xi THEN 1398 
1388 GOTO 1988 

1398 Xi =X 

1488 IF X < 12 THEN 1888 
1418 Ri = Ri + J 

1415 REM == COMPUTE ANGLE = 
1428 IF Al = A3 THEN 1478 
1438 2 = .5 ¥ (Al - A3) /A2 
1448 21 = 14 1/ (2 # 2) 
145 T= -2 * (1 + SOR(Z1)) 
1448 GOTO !4g@ 

1478 T = J 

1488 C=1/ SORC 1 +7 #7) 
1499 S=C#T 

1588 S2=S #5 
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1519 C2=C¥C 
1528 ACK,L) = 8 

1525 REM == TRANSFORM DIAGONAL ELEMENTS == 
1530 AB = 24 AZ *C#S 

1548 ACK,K) = Al #€2 + AQ + 43 ¥ S2 
1550 A(L,L) = Al *# $2 - AB + A3 * C2 
1555 REM == TRANSFORM OFF DIAGONAL ELEMENTS == 
1568 FORT = 170N 

1578 IF 1 < K THEN 1603 

1588 IF 1) K THEN 1648 

1598 GOTO 1748 

1688 AB = ACI,K) 

1618 ACI,K) = C ¥ AB +S ¥ ACI,L) 
1628 ACI,L) = -S # AB + C # ACTL) 
1638 GOTO 1748 

1648 1F 1 < L THEN 1678 

1650 IF 1) L THEN 1718 

1468 GOTO 1748 

1678 AB = ACK,1) 

1688 ACK,]) = € # AQ +S # ACILL) 
1698 ACI,L) = -5 #A@ + C ¥ ACILL) 
1788 GOTO 1748 

1718 A@ = ACK,I) 

1728 ACK,1) = C AB +S # AULT) 
1738 ACL,I) = -$ AB + C # ACL, 
1748 NEXT I 

1745 REN == TRANSFORM MATRIX U TO GENERATE EIGENVECTORS = 
1758 FOR I = 1 TON 

1768 UG = UCI,K) 

1778 UI,K) = C * Ub + S * UCI,L) 
1788 UCL,L) = -S # UB + C # UCI,L) 
1798 NEXT | 

1888 NEXT L 

1818 NEXT K 

1828 RETURN 

1825 REM == NORMALIZE EIGENVECTORS == 
1838 FOR T= 1 70N 

1848 FOR J=1T70N 

1850 U<I,J) = U<I,J) ¥ BCT) 

1868 NEXT J 

1878 NEXT 1 

1888 FOR T= 1T0N 

1898 BCI) = ACI,1) 

1988 NEXT I 

1985 REM == ORDER EIGENSOLUTICN == 
1918 FOR I = 1 TON! 


1928 Jia] +1 

1938 2 = Bil) 

1948 4 = | 

1958 FOR J = 1! TON 


1968 IF 2 ¢ Bid) THEN 1998 
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1978 2 = Bid) 
1988 H= J 

1998 NEXT J 

2888 B(M) = BK) 

2818 B(1) = 2 

2028 FOR J = 1 TON 

2038 2 = Wd,1) 

2048 UCJ,1) = UI,M) 

2058 Wid,M) = 2 

2840 NEXT J 

2078 NEXT | 

2888 RETURN 

3868 PRINT * END OF PROGRAM * 
3818 END 
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SPPEMDIN 8: PROGRAMS USED TO CONVERT VELOCITY HISTORIES 


TE SDeEbEEG TOM Unt S hie 


{8 f mm mmwn emwnn wom mewn em we ww ewww wn wn www ewww sc woe mm www meme wow wn wow wooo on ononn 
28° = PROGRAM #% SIMPS! 4% USED FOR CASE | 

38° INTEGRATES VELOCITY HISTORY WRT TIME TO OBTAIN DEFLECTIONS 

GQ! nnn nn nn 9 on nn no = a on nn nnn nn nnn nnn nn nnn nnn nn nnn oe = oe = 2 oe 
38 CLS : DIM A$<3) : 

68 FC1=8 : FC2)=8 : FC3=8 : FORMIS=" #.RHHH 0°: FORMZ2$=" SHAHHHHA 

70 INPUT" ENTER THE NAME UF ARNOT FIRE ieee seer aeeee gh ULClis cern nn 


88 INPUT’ ENTER THE ORDINATES TO INTEGRATE QVER....cseee 12 "NUMBER : PRINT 
98 INPUT’ ENTER THE SIZE OF THE INTEGRATION INTERVAL...... ",H 
168 CLS 


118 OPEN FILEINS FOR INPUT AS #1; GPEN "DISPR2.CUT* FOR QUTPUT AS #2 


158 INPUT #1,M,TIMECB) VELI1¢@) VEL!2¢8) VEL13¢@) VELSI¢8) VEL62¢@) ,VELS3¢8) 
168 FOR 1 = 1702 

178 INPUT 81M, TIMECI) VELI1C1) ,VELI2¢1) VELI9¢1) VELGIC1) VELS2¢1) VELS3C1) 
188 NEXT | 

198 VIA = VELSI(B) - YELIN(@3 
208 VIB = VELSI{i) - VELII(1) 
218 VIC = VEL6I(2) - VELI1(2) 
228 G1 = (VIA + 4¥WIB + VIC) 
238 V2A = VEL62¢8) - VEL12(@) 
248 V2B = VEL62(1) - VELI241) 
258 Y2C = VELS2¢2) - VEL12¢2) 
260 02 = (V2A + 4¥2B + V2C) #H/ 3 

278 V3A = VEL63¢8) - VEL!3¢3) 

288 V3B = VEL63¢1) - VELI3(1) 

298 V3C = VELS3{2) - YEL13¢2) 

308 G3 = (V3A + 4eV3B + Y3C) #H/ 3 

310 FCI) = Ft) + OF 

328 F(2) = F(2) + Q2 

338 F(3) = FC3) + 93 

348 7 = TIMECI) * 1000 

358 FOR 1 = 1703 

368 IF ABS(QMAX(I)) ) ABS(F(1)) THEN GOTO 380 

378 OMAX(I) = FCI) : THAK(I) = 7 

388 NEXT | 

398 PRINT 42, "T= ‘s: PRINT 82, USING FORMI$:7; 

408 PRINT #2, "nsec DX = ": + PRINT 82, USING FORMZ$;F( 1): 

418 PRINT #2, "in. DY = "ss; PRINT #2, USING FORM2$;F(2): 

428 PRINT #2, "in. 02 = "s: PRINT 82, USING FORM2$:F(3)s:PRINT #2,"in.' 
438 VELII48) = VELI1¢2) + VEL!2¢0) = VELI2(2) + VEL13¢@) = VEL13(2) 


oe 
—< 
ae 
fa 
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440 VEL6168) = VEL6IC2) : VEL6248) = VELS242) : VEL63¢8) = VEL434 2) 
458 IF 4 + 2 ¢ NUMBER THEN SOTO 143 

460 ASCE) = EX) = " 3 A¥C2) = "CY) =": A803) = PZ) = 

478 FORT = 1 10 3 

488 PRINT #2, "THAX"sA$C1); + PRINT 42, USING FORMIS ;TMAX(1): 

489 PRINT 42," msec ODMAX"SASC]) 3: PRINT #2, USING FORN23;GHAX(I): 
4864 RUIN Zeolite 

498 NEXT | 

088 CLOSE #1 : CLOSE #2 


if Sete ee ea eae eee ee ce el ee re ee ee ew ow 
28‘ PROGRAM *# SIMPS2 ## USED FOR CASES I] AND III 

38’ INTEGRATES VELOCITY HISTORY WRT TIME TO OBTAIN DEFLECTIONS 

GQ fo -a-n enna nanan nnn nnn nn nnn nn nnn nnn nnn nnn nnn nnn cnn n ence cnes 
58 CLS : DIM A$¢(3) 

48 FC1)=8 : F(2)28 : FORMIG=" 4H.HNMR oo" : FORMZ$=" HA HAMAAH 


76 INPUT® ENTER THE NAME OF AN INPUT FILE..... ierciaveters sre » ",PILEINS + PRINT 
88 INPUT® ENTER THE GRDINATES TO INTEGRATE OVER... ....een, " NUMBER : PRINT 
96 INPUT" ENTER THE SIZE OF THE INTEGRATION INTERVAL...... au 

188 CLS 


118 OPEN FILEINS FOR INPUT AS #1 : OPEN "DISPR2.OUT* FOR OUTPUT AS #2 
115 PRINT #2,*H = "sHs® SEC 


150 INPUT #1,M,TIMECA) VEL! 208) VEL92(8) VELS2¢8) ,VELI12¢8) VELI92(A) ,VELIS2¢8) 
168 FOR 1 = 1702 

178s INPUT -H1JM, TIMECI) VELI2¢1) VEL92¢1) VELS2¢1) VELII2¢1) VELI92(1) VELIS201) 
{88 NEXT | 

198 VIA = <VELI2¢8) + VEL9208))/2 - VEL52¢8) 

208 VIB = (VELI2¢1) + VEL92(1))/2 - VEL52¢1) 

218 ViC = (VEL12¢2) + YEL92(2))/2 - VEL52¢2) 

228 G1 = (VIA + dB + VIC) #H/ 3 

238 V2A = (VELI12¢8) + VEL192(0))/2 - VELI52¢0) 

248 V2B = (VELI12¢1) + VEL192¢1))/2 - VELIS2¢1) 

258 V2C = (VELI12¢2) + VEL192¢2))/2 - VEL152(2) 

260 02 = (V2A + 4x#J2B + Y2C) ¥#H/ 3 

310 FCI) = FUL) + OF 

328 F(2) = F(2) + 02 

340 T = TIMEC1) * 1008 

358 FOR I = 1 T0 2 

368 IF ABSCOMAX(1)) ) ABS<FLI)) THEN GOTO 390 

378 QMAX(I) = FCI) + TMAX(]) = T 

388 NEXT | 

398 PRINT #2, "T= "s+: PRINT #2, USING FORMI$:7: 

408 PRINT 42, "msec 0¢5,Y) = "; + PRINT 42, USING FORM2$:F( 1); 


9h 


419 PRINT #2, "in.  D(15,Y) = "3: PRINT 42, USING FORM2$:F(2) ; 

428 PRINT 42, ‘in.’ 

436 VEL12(8) = VEL12¢2) + VEL92¢4) = VEL92¢2) : VELS2¢8) = VELS2¢2) 
448 VELI12¢8) = YELI12¢2) : VELI92¢@) = VELI92¢2) + VELI52¢0) = VELi52(2) 
45 IF M+ 3 < NUMBER THEN GOTO 148 

468 AS(1) = "C5,Y) = * : AS(2) = "C15,7) = * 

478 FOR 1 = 1 70 2 

488 = PRINT 42, "THAX":AS(I); + PRINT 42, USING FORMIS ;TMAX(I): 

485 PRINT 42," msec DMAX™:A$(1)+: PRINT 42, USING FORM2$:QMAX(1): 
486 PRINT 42," in.' 

498 NEXT I 

588 CLOSE #1 : CLOSE 42 


SPPENOIN C - TYPICAL CASE 1 ANALYSIS 


in Ehapter VY, the first step 


iL 
a! 
wu 
« 

il 
nt 


Was the construction of an appropriate SAPIW model af the 
equipment substructure under cansideration. ith ls Case, 
the model represents 3 1,808 1b valve attached ta the tree 
end af a cantilever taundatian. 
A. SAPIY INPUT DATA FOR CASE I 

The fallowing data file wat used to generate the moda 
Information tor the Case I model (Ret. 18]. 


HEADING CARD 


CASE 1 -- VALVE ON CANTILEVER FOUNDATION 
MASTER CONTROL CARD 


erg ey ag beg iy Oe, E 


NODAL PT DATA 


Uecker toate 1 hy lk el ang At Se eet ee ala 

Peon 1-6. «94.5.1, 00 

A | er dG. i. 25. 1.8.8 

ce ee |. i ate yas he 4a. 8 
iijeeee —- SD BEAR ELEMENTS 

Pes ie 


MATERIAL PROPERTIES 


ot Me Se aoa. e744 Oly 


a 


Seeman) PROPERTY CARD 


MJ 


1,17.54,9.0,9.8,539.95,32: 


Gi 
NO 
un 

ty 
D 
GS 


tod 


8.0,9.9,8.9,4.8 
ptt gle Clg ta lg tats 
BE oy ae ele saan 
G.8,0.9,9.89,8.8 
BEAM DATA CARDS 
1112.7 01110 :6y8 12 106 


LOMCENTRATED MASS Ob FREE EMD 


Ae oe aaa, Mi ae CUM ge Pie tn Ct eyes 


CYNAMIC AMAL YSIS CARC CIN ACCORDANCE WITH REF.4 


G.09,8, ,0.0,9,8.9 
1,4 
{ft eee on 


EnG OF ole COUMPy TITLE CARD: 
Be 4 St see 


The cutput from the SAPIW cade is nat included here But 
al) the essential information 12 included in the USLUOR code 
eee OCs. 

RB. PICRUST INPUT DATA FOR CASE I! 


The PICRUST code is used to pracecs the SAFIN output 


7 


and calculate influer coefficients carrespanding ta the 


base or support mations of the substructure. Additionall:, 


fu 


cammectivity between the hull and foundation is specified 
ang a substructure file is SGredudced: 


The tollamwing itnput file was used far CASE I: 


Lisi: (OPT Pits 


GM) pe 
ae 
hh 
> 
ba pee 
—s 
— 
ae 
eh pe 
eo 
BQ wi 


MUDAL SPECIFICATION 


3,4 .8 
1 


MODES FROM SAPIY INCLUGED IN ANALYSIS 


EG 
Peat ttitel totter tram time PILRUST code is nat 


imc luded here but all essential cautput from these 


m 


wel atians is included in the USLOE output. 
C. USLOB INPUT DATS FOR CASE I 61414" STANDOFF) 

The folleawing input data wae Bsupplied ta the USLOEB 
code tea specify the parameters for the time integration 


and associated aqutput. A tatal aft L2ul time eteps 


iT! 


oy 


eer ried for the integration interval, each covering an 
Interval af S.GE-4S seconds. The KX*xX parameters denate 
parameters for eq. S which are ammitted due ta the 
eemeat fivity at this equation. 

Tie relating Itput Gata corresccand to the USLOB cutput 


imcluded in this apoendix. 


TIME STEP AND CONTROL CARD FOR INTEGRATION 
ae) oe 14 
EQUATION £5) PARAMETER SPECIFICATION 


3 a BE-& SAK q } S45 q i a “ 
SOUR. mK A yA KOM KKK 


“Cl 
Cn 


BLT REG Gey Weta ls 


G,1,8,8.8,4,8,8,1,.8 
6,9,.1.1,.9,0,8,1,8,8 


SHELL VELOCI Ty SUT PUTS (eat sess Se we 


1,5,14 


SUBSTRUCTURE VELQEI Ty Sielie ites a er Sete eee amie 


a 

Lees peg l «S45 Oa Uy Sawer 

CTRis matation specities the x,y,z tra 
base and tip modes af the subestructur 


chD 


james at the 


ect Ge 


D. LISTING - USLOB QUTPUT FOR CASE ! 
CASE 1 -- VALVE ON FOUNDATION YERT DISPL 


CONTROL PARAMETERS FOR SUBSTRUCTURE =| 


NPOINT = 7 NDOF = 18 
NCOOF = 3 NUDOF = = 15 
NFREG = = 18 NFACE = i 

NSYMTY = 8 NFORCE = = 13 


CONNECTION DATA FOR INTERFACE POINTS 


SHL = SUR x it Cee Wl eece Coorl LEN) DEGCON 
36 i 1 1 1 8 8 8 3018 27 2.78088RE+E2 
GDEGROT = 9. B8888RE+81 
JEQUATICN NUMBERS FOR SUBSTRUCTURE 1 
POINT xX Y 2 XK W 22 
feeieee co Ser 8B 8 
Peeteeo- eo te BCC 
Seo - 9 ete 8 8 
4 48 wf 12 8 @ 8 
Sicld 14 13. 8 8 28 
Stites i 0 6 
ete oe 70 ee BC 


FIXED-BASE NATURAL FREQUENCIES FOR SUBSTRUCTURE |! 


FREG 
1.34387E+83 
1. 845495404 


MODE FREQ 
4 4,93785E483 
8 1.21696E+84 


MODE FREQ 
ASS41E +82 
68871E+83 


.22848E+84 


MODE FREQ 
2 1,85485E+83 
6 9.55885E+83 


18 1.88142E+84 


MODE 
d 


= 


1 4 
a ai? / 
yet 


ICOMPACTED FIXED-BASE MODES FOR SUBSTRUCTURE 1 


| 


2 3 


4 i 


Sy puceloes) .coteeoet 
ALS6BE-28 8.97 11E-89 
2402E-81 -4.8588E-28 
»BI8S6E-18 -2.4745E-81 
.2987E-28 = 1./988E-88 
4745E-81 -9.6785E-28 
pPBOZESt Gera. G7) Lend 


et 
eee 
Sia! 
43 
3 4 
Onze 
= 
8 9.4185E-28 2.6743E-08 


Pease oOmrc,. 2ovverle 82.47 / dc) 
PeC4OcESUlge liesceic 1 Of arcs 
Bo /LS7Ecil wl l47i7etOe =2,.0/05Ete0 
POSSE SUe SIR OSCE Sl a. 7chces! | 
2,4745E-0)) },ouboal2 2./dede7i3 
J697B2E-11 -2,3685E+88 -1.4497E+88 
2,8465E-89 9.788BE-12 3.9373E-11 
Geos Mesti@elte/OlLG-12 wevolicEatS 


9? 


9 -Z.6971E-G1 -].4468E=19 S.e7S7ba taco sy cette ee oe anne 
16 7.5646E-18 -4.9020E-81 3.7793E709 5707 sin 2 a2 cote 
1] {.2878E-17 S.S9S7E-08° S228 28E=01 “SENG eo Zee cwe seas 
12 4, 9B28E-81 <1 PI73E-12 3.5) 72E-s) aime cette tov cetre 
13. 9. S88BE-19 -6,8836E-01 4. 6B405-07 —soeodseai gmetad ice 
{4 1,9485E-19 4.4806E-88 6.8836E-8! -3.9949E-14 -1.8775E-14 
15 =6.8636E-81 <2.3799E-17 —2, 3536E aie 9) Goer te see eects 

é 7 8 9 18 
LL .4749E+88 3.5428E-14 1.2233E-28 1.2895E-12 -2.3705E+88 
2 1.9626E-13 2.2993E-11 -1.4719E+08 6.6505E-13 9.7987E-11 
3 2,4783E-18 -2.3685E+88 -3.3464E-12 -1.4634E+88 -1.8812E-12 
4 2,36856+88 S.384BE-14 1.5886E-28 -7,d508E-13 -1.4497E+88 
2 28925E-13 3.697711 --2.3685Et00 VinececEo 2 oe coe 
6 ~{,0816E-18 1.47196488 S.45B1E-12 2.36936+88 -3.7449E-12 
7 =2,3393E+88 4.9358E-14 1.1252E-28 -4.6836E-12 1.4839E+80 
8B 2.1650E=13 256543E—1) <2 ea0sb BO letaloe 12 senscace a 
9 -1.4869E-18 1.4538E488 -5.4572E-12 -2.37246+88 -3.6047E-12 
18 «1 .3757E+08 22 7551ES14 4.47 1GEo 2 aera Vor oes 2esogee as 
{i {,2928E=13 2. 188SE=1 <f.3757Et80 See2a77E-12 eaves 
13 2.4276E-19 -2.3754E+88 3.3858E-12 1.4715E+88 =7 742g9E te 
16-9, 3383E-82 -1.9932E-19 -4.79308-22 =, 2783E—12 da edoce awe 
13 -8.6446E-15 -1.43988E-12 9. 3383E-§2 =4,1796E-14 =7776ie oe 
16 =2.291GE=12 2.24286-02 <2. 3083R=1a) =| ObgoGa2 =< eleccee 


{SUBSTRUCTURE CONTROL PARAMETERS 


SU8 NFRE@ NDOF NCDOF NPTS LOC@ NEG NBLOCK HROWFC 
| 12 18 3 i wale 13 i 3 
N8LFOR KOLBAR N8LBAR SEGFR@ 8EGSSL 8EGBAR S8EGFOR 
$254 | | | | i 
CONTROL PARAMETERS FOR SOLUTION 
NTIME = 1281 NEQS = 638 
NSKIP = 2 NRECS = 681 
NCHRG = | NTHETA = =|? 
NQUAD = 9 NDELTA = = 21 
NFINE= 1 N@D07 = 264 
KOUPLE = 2 NS5U85 = | 
SOLUTION PARAMETERS 
TOTALM = 3,541492E+04 AREA = 3.639844E+86 
RHOFLU = 9.452168E-85 CFLU = 6. 800888E+84 
DELT = 3.888888E-86 TEND = 6,088008E-83 
ALOAD = 8.844888E+82 STOFF = 1.4/4800E+83 
RLOAD = 1,686888E+83 ORNRAD = 8, 


98 


DATA FCR EMPIRICAL PRESSURE-TIME HISTORY 


WEHRG = 
PeMLT = 
PZEXP = 


2. d88008E+83 
KXAAXAXAKAA 
XXEXAKAXKAKK 


THEXP = 


IVELOCITY DISTRIBUTION AT SHELL STATION 


Ser 


TIME 


8. 

1. 0888E-85 
2, 8888E-85 
3. 8888E-89 
4, 8808E-85 
3. B8B8E-85 
6,0888E-85 
7. 8888E-89 
8.8886E-85 
9, B888E-85 
1. 8088E-84 
1. 1888E-84 
. 2808E-84 
» 3888E-B4 
4888E-84 
 JOBBE-84 
 6888E-84 
./8B8E-84 
. S8BBE-84 
 98BRE-84 
2, 0888E-84 
2.1888E-84 
2, 2888E-84 
2, 3808E-284 
2.4800E-84 
2, 9868E-84 
2.4888E-84 
2.7B888E-84 
2,3806E-64 
2, 7808E-84 
3. 8808E-84 
3.1886E-84 
3. 2880E-84 
3. 3888E-84 
3, 40B8E-84 
3. 3808E-84 
3,6888E-84 
3. 7B88E-84 
3. 8808E-84 


ome Pema ei eee fee femelle 


WOOT! 8) 


8, 
4. 94446-81 
1. 8855E+868 
3.8391E+88 
4,1484E+88 
8, 4881E+88 
1.1347E+81 
1.41456+81 
1,7883E+81 
2.0102E+81 
2. 3241E+81 
2.4947E+81 
2.9933E+81 
3, 3481E+81 
3.7123E+81 
4,8922E+81 
4,4771E+81 
4.8745E+81 
Se2/2leeet 
5.4819E+81 
6, A891E+81 
6,5836E+81 
6, 91616481 
7. 3292E+81 
7.7437E+81 
8,1493E+61 
8,54606E+61 
8,9472E+91 
9, 34469E+81 
9,7634E+81 
1] ,8154E+82 
1,0542E+82 
1,0933E+82 
1,1317E+82 
1,1781£+82 
1. 2886+82 
1.24465E+82 
1,2844E+82 
1.3226E+82 


WOOT. 98) 


8, 

7. 0446E-82 
2 .8722E-81 
3.8428E-81 
9,3465E-81 
1. 3193E+88 
1,7214E+88 
2.14166+88 
2.57 92E +88 
3,8245E+88 
3.4895E+88 
3.9497E+88 
4.45 93E+068 
4, 94516+08 
3.47846E+68 
6, 88545+08 
§,5359E+88 
7 .8787E+88 
7.6848E+88 
8,1394E+68 
8.4490E+88 
9, 1899E+88 
9, 7044E+88 
1,0209E+81 
1.0788E+81 
1.11886+81 
1,1643E+81 
1,2097E+81 
1.2528E+81 
1.2918E+81 
1,3381E+81 
1,346469E+81 
1.4826E+81 
1.4379£+81 
1.4712E+83 
1,5834E+81 
1,.5354E+81 
1,5655E+81 
1,59948E+61 


7 


36 


SURCUT = XXAXAXKAXXAXX 
THALT =  AXAXXAAKXXAX 
KXAAKXAXKKAKA 


(SEGMENT 


WOOT? 188) 


8, 
1.7639E-81 
6.5644E-81 
1.30056+88 
2.8219E+88 
2,7649E+88 
3.4983E+88 
4.2275E+88 
4, 95596 +08 
3.46645E+88 
6.3713E+88 
7 8994E+8B 
7 ,.8237E+88 
8.5456E+868 
9, 3141E+88 
1 8891E+81 
1.8856E+81 
1.1647E+81 
1.241 1€+61 
1. 3288E+81 
1,3957E+81 
1.4723E+81 
1.5459E +81 
1.6185€+81 
1. 6899E+81 
1.7567E+81 
{.8247E+81 
1,8895E+81 
1. 9589E+81 
2.0186E+81 
2.8449E+81 
2.1219E+81 
2.1774E+81 
2.2308E+8 | 
2.2828E+81 
2. 33086481 
2. 30086481 
2.4379E+81 
2.4892E+61 


So PGINT 


WOOT 278) 


8. 
7 445E-82 
2.8728E-2! 
3.84186-81 
9, 3434E-81 
1.31866+88 
1,7282E+88 
2.1397E+88 
2.5743E+88 
3.8224E+88 
3.4840E+88 
3.9625E+88 
4, 4502E+80 
4,9539E+88 
3.4658E+88 
5. 9871E+88 
§.51467E+88 
7 .8482E+88 
7.07 88E+88 
8.1897E+08 
8.4353E+88 
9 ASIFE+88 
.6618E+88 
OL S1E+81 
04486481 
11226481 
1588E+81 
»2027E481 
2451 E+81 
28300481 
»S2O9E+ 81 
do BE+B1 
 S918E+81 
1.4261E+81 
1,4585E+81 
1,4894E+81 
{.5205E+81 
1,5494E+81 
1,9767E+81 


te ema eee enema fete fmt eee tet 


19) 


131 
133 
133 
ay 
133 
16! 
163 
165 
167 
169 
171 
173 
Mo 
177 
179 


3, 9808E-84 
4, 8000E-84 
4, 1880E-84 
4, 2808E-84 
4, 3088E-84 
4,4880E-04 
4. S086E-84 
4, 6800E-84 
4. 7088E-84 
4, d800E-84 
4. 98B0E-84 
3 B888E-84 
3. 1G0RE-84 
3. 200BE-84 
3 3888E-84 
3. 4008E-84 
3. J8B8E-84 
3, 6888E-84 
3./80GE-64 
3 88GGE-84 
3. 7888E-84 
6. B808E-84 
6, 1008E-84 
6, 2008E-84 
6, 3808E-84 
6,4800E-84 
6, 5080E-84 
6, 4808E-84 
6./888E-84 
6, 8800-84 
6, 9888E-84 
7 BG8BE-84 
7 AGGGE-84 
7, 2808E-84 
7. 3888E-84 
7, 4808E-84 
7, S88GE-84 
7, 6880E-84 
7, 7888E-84 
7, S888E-84 
7, 9BQGE-84 
§.8800E-84 
8.1888E-84 
8. 2008E-84 
8. 3808E-84 
§.4888E-84 
8, 9080E-84 
8. 6088E-84 
§.7888E-84 
8. 8080E-84 
8. 9088E-04 


1, 35995482 
1.3974E+82 
.4349E+82 
A719E+B2 
, JO86E+82 
» S454E+82 
O81 7E+82 
»6175E+82 
1.6533E+82 
1,6888E+82 
1. 7237E+82 
1.7588E+82 
1,7928E+82 
1, 8255E+82 
1,8585E+82 
1, 8910E+82 
1.9229E+82 
1, 9548E+82 
1, 98456+82 
2.0147E+82 
2.8442E+82 
2.0738E+82 
2,1818E+82 
2.1 283E+82 
2. 1558E+82 
2,1811E+82 
2,2867E+82 
2,231 7E+82 
2.2063E+82 
2.2885E+82 
2. 3041462 
2.3271 E+82 
2. 3497E+82 
2,371 9E+82 
2.3936E+82 
2.4147E+82 
2 .4354E+82 
2.4558E+82 
2.4757E+82 
2, 4958E+82 
2 D1 39E+82 
2.0320 +82 
2. D0 8E+B2 
2. 0685E +82 
2. 9857E+82 
2.6826E+82 
2.6193E+82 
2,.6350E+82 
2.601 8E+82 
2.6663E+82 
2,6813E+62 


feminist ei tee 


1,6222E+81 
1, 6490E+81 
1.6748E+81 
1,469787E+81 
1.7222E+84 
1.7432E+61 
1.7639E+81 
1.7843E+81 
1,8817E+8 
1, 81806+81 
1.8342E+81 
1.8487E+81 
1, 8681E+81 
1.8714E+81 
1,8824E+61 
1. 8988E+81 
1.8957E+61 
1. 9813E+61 
1,9863E+61 
1. 9671E+81 
1.9874E+81 
1.9876E+81 
1.9857E+81 
1,.8995E+61 
1 .8932E+61 
1, 8869E+61 
1,8798E+61 
1, 8681E+81 
1.8572E+81 
1.8463E+81 
1, 8348E+61 
1,.81896+61 
1, 8835E+61 
1.7882E+681 
1.7723E+81 
1, 7531484 
1.7336E+81 
1.7146E+81 
1.6961E+81 
1.46736E+61 
1,6584E+61 
1.6279E+81 
1.6862E+81 
1, 5889€+81 
1. 9546E+64 
1,.5298E+81 
158396481 
1.4778E+8} 
1.44885 +61 
1,49208E+81 
1,391 7E+81 


168 


2.038ce +81 
2. 2884E+81 
2.6395E+61 
2, 6885E+81 
2,7371E+81 
2.7866E+81 
2 C308E+ 81 
2.0822E+61 
2. 9384E+8! 
2,9778E+81 
3.8233E+61 
3,8677E+61 
3.4125€+81 
3.1558E+81 
3.1964E+81 
3.2371E+61 
3,2765E+81 
3,3136E+81 
3.3489E+61 
3,3847E+61 
3.4182E+81 
3.4496E+81 
3.4789E+81 
3,5074E+81 
3. 0335E+81 
3,5083E+81 
3.5823E+81 
3. 6856E+61 
3.6277E+81 
3.6484E+81 
3.6682E+61 
3,6873E+81 
3,7054E+81 
3./7222E+61 
3.7379E+84 
3. 7534E+81 
3.7676E+61 
3,/7886E+81 
3.7924E+81 
3, 8842E+81 
3.8147E+8! 
3,8242E+81 
3.8326E+8! 
3,8407E+81 
3. 8481E+81 
3.9542E+81 
3.8594E+81 
3,8640E+61 
3,8680E+61 
3.8716E+61 
3.8736E+81 


1, 6836E+8! 
1,62916+01 
1.6526E+81 
1,6759E+83 
1.6979E+81 
1.7174E+83 
1,7365€+81 
1,7553E+81 
1,7718E+84 
1,7856E+61 
1, 8888E+81 
1,8127E+81 
1 .8223E+61 
1.8317E+81 
1, 8468E+81 
1,8463E+81 
1 8581E+81 
1, 8537E+81 
1.8565E+61 
1,8552E+81 
1 ,8533E+ 81 
1,851 2E+81 
1.84716+61 
1,8386E+81 
1, 8308E+81 
1,8213E+81 
1, 81185481 
1,7978E+81 
1 .7844E+81 
1,7711E+04 


1, 7564E+81 


1.7389E+81 
1. 7211E+61 
1.7034E+81 
1. 6851E+81 
1,4635E+61 
1, 6415€+01 
1.6281E+61 
1,5993E+61 
1.5744E+81 
1.5489E+64 
1,5240E+81 
1,5081E+63 
1.4725E +84 
1.4448E+01 
1.4161E+81 
1 .3888E+81 
1,.3686E+81 
1.3294E+84 
1. 2985E+8! 
1,2682E+81 


13} 
183 
185 
187 
139 
19] 
193 
195 
177 
ey 
281 
283 
289 
287 
209 
214 
Z13 
ihe 
217 
219 
221 
223 
(es 
227 
229 
231 
233 
233 
237 
239 
241 
243 
245 
247 
249 
201 
206 
doe 


28! 


9, 8888E-84 
9, 1808E-84 
9, 2098E-84 
7, 3eg8E-84 
9, 4808-84 
9, S888E-84 
9, 6808E-84 
9, 7BOBE-94 
9, BBBGE-94 
9. 9GAGE-84 
1 .6888E-83 
1.8198E-83 
1.9288E-83 
1.8388E-83 
1 .9498E-83 
1,8588E-83 
1.8689E-93 
1.8 7B8E-83 
1. 8888E-83 
1. 8988-83 
1.1888E-83 
1, 1188&-83 
1.1288E-83 
1,1388E-93 
1.1498E-93 
1.1588E-83 
1 .1689E-83 
1.1 788E-83 
1.1888E-83 
1.1998E-83 
1. 2888E-83 
1.21 88E-83 
1.2289E-83 
1, 2388£-83 
1.2499E-83 
1,2588E-83 
1. 2688-93 
1.2798E-33 
1.2899E-93 
1. 2988E-83 
1.38989E-83 
1.31 886-83 
1.3208E-83 
1.3389E-93 
1,34988E-83 
1. 3588E-83 
1.3688E-83 
1.3788E-83 
1, 3388E-83 
1, 3999E-83 
1. 4808E-83 


/OIS1E+O2 
7182E+82 
7 238E+82 
7372E+82 
» /O83E+92 
7633E+82 
feoge+O2 
.7873E+82 
AIFIE+G2 
2.8185E+82 
2,821 6E+82 
2.8328E+82 
2 .8423E+82 
2.8524E+82 
2 .8623E+82 
2.8715E+82 
2 ,8892E+82 
2.8887E+82 
2 .8971E+82 
2. 9854E+ 82 
2.91 38E+82 
2.,9281E+82 
2.9272E+82 
2.9342E+82 
2.941 1E+82 
2.9473E+92 
2. 9531E+82 
2. 9589E+82 
2.9646E+82 
2.9793E+82 
2.9752E+82 
2.9796E+82 
2. 9841E+82 
2. 9886E+92 
2.9929E+82 
2.9966E+82 
2.9999E4+92 
3.9831E+82 
3.8863E+82 
3, 8094E+82 
3.9121£+82 


z 
2 
2 
2 
2 
2 
2 
2 
2 


1.3643E+81 
1,3335E+8! 
1,3019E+8! 
1.2789E+81 
1.2484E+81 
1,.2198E+81 
1.1765€+91 
1,1429E+81 
1,1198E+91 
1,0778E+81 
1.0442E+81 
1, 9889E+91 
9, 7388E+98 
9, 3938E+99 
9, 8548E+98 
8.7844E+88 
8. 2494E+08 
7. 9933E+88 
7.6462E+88 
7. 3836E+88 
6.94] 3E+88 
6.5667E+88 
§.1998E+88 
5. 8414E+88 
5.4887E+88 
J+1148E+88 
4.7374E+88 
4, 3669E+88 
4 ,0038E+80 
3.6475E+98 
3.271 1£+98 
2. 8911E+88 
2.9193E+88 
2.1528E+88 
1.7896E+98 
1.41468E+88 
1 .9388E+88 
6,6715E-81 
2. 9858E-8 1 


sOnasGae-O7 
-4,3242E-84 


3.8142E+82 -8,8686E-9! 
3.8164E+82 -1.1766E+08 
3,9185E+82 -1.5432E+88 
3.8286E+82 -1.9846E+09 


3.8225E+82 


-2,.2646E+88 


3.8236E+82 -2.4328E+08 
3.9246E+82 -2,9962E+99 
3.9256E+82 -3.3579E+88 


3,8267E+82 


-3.7125E+68 


3.0278E+82 -4.8645E+88 


19] 


3,8739E+81 
3.3746E+81 
3.8746E+81 
3.8734E+8! 
3,8713E+8 
3.8686E+8 1 
3,8652E+81 
3,8688E+8 1 
3,.8952E+81 
3.8489E+81 
3,8421E+8! 


3,.3346E+81 © 


3,8257E+81 
3.81463E+8! 
3, 6061E+81 
3.7955E+8! 
3,7849E+81 
3,7734E+81 
3.753 1E+8! 
3.7488E+8! 
3.7345E+8! 
3.7205E+8! 
3.7858E+8! 
3,6898E+81 
3.6729E+81 
3.6955E+8! 
3.6375E+8! 
3.6185E+81 
3.5985E+81 
3,5777E+8! 
3,5971E+8! 
3, J399E+8! 
3.9138E+81 
3,4909E+81 
3.4676E+81 
3.4442E+81 
3.4205E+8! 
3, 3959E+8! 
3,3707E+8! 
3.3451E+8! 
3, 3198E+8 1 
3.2925E+81 
3.2647E+91 
3.2363E+ 81 
3.2070E+8! 
3.1774E+81 
3.1478E+8! 
3.1175E+8! 
3.9869E+8! 
3.9553E+81 
3.9233E+81 


1,23886+81 
1, 2048E+81 
1,1725E+81 
1,1397E+91 
1 JO77E+8! 
1.9754E+81 
1,0493E+61 
1,6952E+81 
9. 7897E+88 
9,.3739E+88 
9 ,8253E+88 
3,6604E+98 
8. 2991£+98 
7. 9449E+98 
7 95 7E+88 
7, 2292E+88 
6, 87496498 
§,3125E+88 
6,1599E+88 
3,81 28E+88 
3. 4468E+88 
5. 8596E+88 
4 .7009E+98 
4,3419E+88 
3,9896E+88 
3.6193E+98 
3.2426E+88 
2.8752E+88 
2.91 83E+88 
2,1684E+88 
1,7998E+88 
1.4288E+88 
1.8673E+068 
7.4232E-91 
3.6192E-81 
3,2835E-83 
-3,5943E-8! 
~7 1455E-81 
~| ,0653E+88 
-1,.4187E+88 
-1.7571E+88 
~2, LIB1E+88 
-2,4572E+98 
-2, 8808E+88 
~3,1364E+99 
-3,.4782E+99 
-3,8184E+88 
~4,1462E+88 
~4,4775E+98 
-4,8823E+88 
-5.1225E+88 


283 
289 
287 
289 
291 
293 
fase 
297 
299 
381 


41 RBE-83 
4208E-83 
4308E-83 
44885-83 
 4508E-83 
4688E-33 
4788E-83 
4888E-93 
4988E-83 
 J088E-83 


eerie eee eee et 


3.8288E+82 
3,8238E+82 
3.8279E+82 
3,8288E+82 
3.8281E+82 
3.8277E+82 
3.8267E+82 
3,8257E+82 
a 


-4,4214E+08 
-4,7756E+88 
-5.1262E+88 
-3.4737E+88 
-3.8188E+88 
-S, 1611E+88 
-$,9018E+88 
-§.8489E+88 


O247E+82 -7,1776E+88 


3,8239E+82 


-7 91216 +88 


2,9917E+81 
2. 9599F 38} 
2.9276E+81 
2.8944E+81 
2,0606E +81 
2.8268E+81 
2,7936E+81 
2,7 099E+8} 
2.7259E +8} 
2,6903E+91 


-3.4443E+88 
-3.7542E+88 
-§ 081 3E+98 
-3,3720E+68 
-$,6982E+88 
-7 ,8821E+88 
-7 ,Jd24E+ 88 
-7 ,9999E+88 
7, O941E+08 
“8, 1858E+83 


( This history is abbreviated for illustrative purposes.) 


WVELCCITY DISTRIBUTION ALONG SU8STRUCTURE 


STEP 


“wt of wo 


TIHE 
Q 


VELG 
8, 


1,1) VEL 


1,8888E-85 -7.46928E-84 -6.4883E-83 


2.8888E-85 -3.6167E-83 -2.3780E-82 
3.8888E-85 -8.8848E-83 -3,3246E-82 
4. OG8BE-85 -1.6429E-82 -3.7845E-82 
D.8888E-85 -2,6541E-82 -2,6866E-82 


6.8808E-85 
7, 8888E-85 
8. 8888E-85 
9, 8888E-85 
1, 8888E-84 
1. 1888E-84 
1. 2888E-84 
1, 3008E-84 
1.4888E-84 
1. 5888E-84 
1, 6888E-84 
1.7888E-84 
1 ,8888E-84 
1, 9888E-84 
2. 0888E-84 
2. 1888E-84 
2, 2808E-84 
2, 3888E-84 
2.4880E-84 
2. 0888E-84 
2.6808E-84 
2./888E-84 
2.8808E-84 
2, 7888E-94 
3, 8888E-84 
3,1888E-84 
3.2088E-84 
3,3888E-84 
3.4888E-84 


-3.8728E-82 


-5,2807E-82 


-$,8819E-82 
“8, d171E-82 
Sie Sg2ZE—ni 
-1,2514E-81 
-1,4552E-81 
-1,6745E-81 
+), 02Z7E 6H 
S2igeleedn 
-2,3783E-81 
-2.46285E-6! 
-2,8685E-81 
“3, 1149E-81 
-3,3681E-81 
=. oe 7e<08 
-3, 8688E-31 
~4,1029E-81 
~4,3351E-81 
-4.9982E-81 
-4,7767E-81 
-4.9885E-81 
= Onto aeeadll 
-3,d628E-81 
SP cess 
3. 7835E-01 
-3,8558E-81 
-3.9947E-81 
-$.1224E-81 


3.8282E-84 
4, 4051E-82 
1.4789E-81 
1,8835E-81 
2.9143E-81 
4, 1250E-81 
J. Ji 22E-81 
7 1056E-81 
8,8989E-81 
1,0795E+88 
1.2921E+88 
1,5173E+80 
1.7626E+80 
2.0222E+88 
2.2978E+88 
2, 9879E+ 88 
2. 8949E+ 88 
3,.2141E+88 
3. 59481E+88 
3.8942E+08 
4,2558E+88 
4,6265E+08 
J. 8121E+88 
3.4117E+08 
5.8235E+08 
6, 2440E+00 
6,46781E+83 
7 LI91E+88 
7, 9694E+88 


182 


1 


VEL¢ 
8, 


2,8/728E-8! 
2 8418E-8! 
9,3436E-81 
1,3186E+88 
1.7282E+88 
2.1397E+88 
2.0763E+88 
3.8224E+88 
3.48485 +88 
3.9625E+88 
4.45825 +88 
4,95 39E+ 88 
2. 4698E+88 
2. 9891E+88 
6,5167E+88 
7 .8482E+88 
7.9788E+83 
8.1897E+88 
8.6353E+88 
9, LSL9E+88 
9,66:8E+38 
1. 8161E+81 
1,8648E+81 
Pei lZce+ ei 
1,1588E+81 
1.2827E+81 
1,2451E+8] 
1.2835E+8! 
1.3287E+81 
1,3578E+81 
1.3918E+8) 
1.4261E+81 
1,4585E+81 


VEL¢ 
Q 


b,1) 


7. O445E-82 §.6872E-07 


1. 9828E-86 

1. d989E-864 
=eaglwekcto 
~4,3943E-86 
-2,8589E-85 
-1,.2167E-84 
-3,7888E-84 
-8, 862HE-84 
-1,7841E-83 
-3.1668E-82 
=. le=es 
(seus 
1,181 7E-82 
-1,9895E-82 
=l ye 72 eeauc 
“200 0ETOL 
-J,.2914E-82 
4 aac yeaee 
4, 7349E-82 
-5,9691E-82 
=Paecie-O2 
-8 .4681E-82 
“9, 9248E-82 
-1,1536E-81 
Toe vESe! 
Sede cman 
=), 7266E=81 
=Ta7seaelt 
Seq Ge7 cate 
-2,4339E-81 
-2.7082E-81 
-2,981 2E-81 
See? Cue en 


VEL( 
0. 
=1,7765E-89 
23,07 74E700 
1.2365E-69 
4. 4999E-85 
Sine 7 es04 
-9,1328E-84 
-2,3448E-83 
-4,2651E-83 
-§ .8883E-83 
Sonori deste 
SO, dr teste 
=le2dScboue 
7 .8835E-83 
2.8128E-82 
3,8738E-82 
Org) secu 
Piey Eau 
1.4836E-81 
peayzeool 
2.6265E-81 
3.4499E-81 
4,4315E-81 
9. J824E-81 
6,9154E-81 
8, 4451E-81 
1,0198E+68 
1.2161E+88 
1.4379E+868 
1.6854E+88 
1. 9573E+88 
2.2683E+88 
2, 9898E+88 
2.7461E+88 
3.3327E+88 


é,2) 


VEL¢ 
a. 
-4,3588E-14 
1.4181E-12 
S.e1s0bsie 
1. 3085E-88 
1.8281E-87 
1. 3878E-86 
8.6295E-86 
3.7798E-85 


5,3) 


1.8797E-82 
3.1738E-82 
3+84528E-82 
7,.6605E-82 
1.1865E-81 
1 .53d7E sa 
2,8506E-81 
2.607 3E-8) 
Sioce/ tour 
4.4397E-81 
9, 8061E-81 
dev evckaey 
7,81 95E-81 
8.1594E-81 
9 3780-81 
.0733E+88 
2163E+88 
. 3685E +88 
JLIFE+ 8B 
. 7BB1E+88 


eee het 


147 
149 
13] 
Io 
139 
13? 
[oY 
16] 
163 
165 
167 
169 
17! 


3, 3888E-845 
3.6888E-84 
3, 7BBBE-84 
3. 8088E-84 
3.9888E-84 
4, 8888E-84 
4, 1888-84 
4, 2888E-84 
4,38806-84 
4.4000E-84 
4, 5888E-84 
4. 6888E-84 
4, 7888E-84 
4, 8888E-84 
4, 9989E-84 
2. 8888E-84 
2.1896E-64 
2. 2088E-84 
3. 3886E-84 
3. 4888E-84 
2. J8B8E-84 
2. 6888E-84 
J. 7888E-84 
J. 8888E-84 
2. 98BBE-84 
6.8808E-84 
6, 1888E-84 
6. 2888E-84 
6. 3888E-284 
6.4888E-84 
6. 888E-84 
6.6888E-84 
6./888E-84 
6, 8888E-84 
6, 9888E-84 
7 BBBBE-84 
7, 1888-84 
7. 2808E-84 
7, 3888E-84 
7. 4888E-84 
7 S888E-84 
7. 6888E-84 
7, 7A88E-84 
7 .BOG8E-84 


-5, 2468E-8! 
-6,. 3318-8! 
-6,4496E-61 
~$.5416E-81 
eoroe27E—bi 
-6.7834E-81 
"6. 77396E-81 
-4,8457E-81 
6. 9893E-81 
-6.9747E-81 
-7,8376E-81 
-7, 1828E-81 
-7.1622E-81 
my ecole tl 
Secovlc =O) 
e(ado17E>0! 
7 .4298E-81 
-7 .5894E-81 
-7, 5834-81 
-7,6648E-81 
Sve doorce! 
Sy Colle?e! 
aac oc- 8! 
-8.8498E-81 
-8.1562E-81 
-8,2643E-81 
“8. 3884E-81 
-8.5826E-81 
=enocbde 70) 
-8.7428E-81 
-8,.87B1E-81 
B77 00E-81 
ee etecoaeO! 
piace sl7E-0) 
arma Y7eou 
~9.5874E-81 
-9 6374-81 
aye 7ece-Ol 
avec) dee-Bl 
-1,B821E+88 
-1,B158E+88 
-1.8277E+88 
-1,3396E+88 
-1.6519E+88 


7. 98BB8E-84 -1,8647E+88 


8.8808E-84 


-1.8765E+88 


8.1888E-84 -1.8881E+88 


8. 2888E-84 


-1.1882E+88 


8. 3888E-84 -1.1121E+88 


8.4888E-84 


SZ JeetO0 


8.5888E-84 -1,.1348E+88 


3,02006+98 
8, 4948E+88 
8.9486E+89 
9, 4509E+88 
9, 9375E+88 
1. 8431£+81 
1, B931E+81 
1.1437E+81 
1.1948E+81 
1.24465€+81 
1. 2989E+81 
1901 2E+81 
.4043E+81 
4588E+61 
SIISE+81 
SS59E+B1 
.6208E+81 
.6768E+81 
1.7311E+6! 
1.7871£+81 
1.8433E+81 
1,8995E+81 
1. 9559E+81 
2.8129E+81 
2.8782E+81 
2.1275E+81 
2. 1851E+83 
2. 2436E+81 
2.3023E+81 
2. 3608E+81 
2.4194E+81 
2.4788E+81 
2, 0388E+81 
2,097 1E+81 
2.606248 
2.7162E+81 
2.7768E+81 
2.8354E+81 
2.8958E+81 
2.9502E+81 
3.8155E+81 
3.8755E+81 
3.1352E+81 
3.1957E+81 
3.2563E+81 
3.3167E+81 
3.3766E+81 
3.4372E+81 
3.498BE+81 
3,9086E+81 
3.6189E+81 


eee et 


183 


1.4894E+8! 
1.5285E+61 
i O494E+a} 
1.5767E+81 
1, 6836E+81 
1.4291E+81 
1,46526E+81 
1.6759E+81 
1,.6979E+81 
1.7174E+81 
1,7365E+81 
1.7593E+8! 
1,7718E+81 
1.7856E+81 
1. 8808E+81 
1,.8127E+81 
1,8223E+81 
1,8317E+81 
1.8488E+81 
1,.8463E+81 
1. 8581E+81 
1 .8537E+81 
1 ,.8565E+81 
1.8552E+81 
1,8533E+81 
1.8512E+81 
1,8471E+81 
1 .8386E+81 
1. 8388E+81 
1.8213E+81 
1.d118E+8! 
1.7978E+8! 
1.78445 +81 
1. 7711E+81 
1.7564E+81 
{,7389E+81 
Fey ciete) 
1. 7034E+81 
1, 6851E+81 
1.46635E+81 
1, 64156+81 
1.6281E+81 
1.5993E+81 
1.5744E+81 
1.5489E+81 
1,.5248E+81 
1. S0B1E+81 
1.4725E+81 
1.4448E+81 
1.4161E+81 
1, 3888E+81 


-3.9835E-6! 
-3,7826E-81 
aecdesee) 
S407 Loe= 8) 
-4, 9185E-81 
sou? dorcel 
-5.347E-31 
-6, 8813-81 
oma eae a0 
-6,/7464E-81 
aval ely e—al 
-7.4965E-61 
-7 B69 1E-8! 
8, 2381£-81 
-8.6922E-81 
-8.9682E-81 
“9, 3189E-81 
Sica Job =0 
Senor Lest. 
1 O311E+88 
-1,.0623E+98 
-1,0923E+88 
-1,1289E+88 
-1.1488E+83 
aely coer ee 
~1.1974E+88 
Spel y7Edue 
-1.2403E+88 
-1,2592E+88 
-1,2764E+88 
-1. 29196 +86 
-1,.3857E+88 
-1,3176E+88 
eieder7ered 
1. 3348E+88 
eladeoetee 
-1,3473E+88 
-1.3584E+88 
1.35196 +88 
-1.3518E+88 
-1.3582E+88 
-1.3472E+88 
1. 3429E+88 
-1.33716+68 
-1, 33818 +88 
-1.3219E+68 
1.31256 +88 
-1.38216+08 
1, 2989E+68 
-1.2789E+88 
-1.2664E +88 


3, 7494E+90 
4,19496+88 
4,5753E+88 
9.1844E+88 
5./7236E+88 
6, 2922E+88 
6,8896E+86 
7, 5155E+86 
8.1493E+88 
8 ,3503E+88 
9,5574E+88 
O2B9E+B 1 
/1B44E+91 
AABISE+OS 
261 4E+81 
»S426E+81 
.4253E+8 1 
2095E+8 | 
DIOBE TBI 
»O814E+8! 
.7686E+81 
»BI62E+B1 
944 2E +81 
2. 8321E+81 
2,1199E+81 
2. 20 74E+8! 
2.2944E+81 
2, 3086481 
2.4659E431 
2 J499E+84 
2, 6320E+6 1 
2.713é6E+81 
2.7938E+81 
2.8784E+83 
2.9464E+81 
3, 62B1E+81 
3.6916E+81 
3.1688E+8! 
3.2277E+81 
3.2922E+81 
3.3543E+81 
3.4148E+81 
3.4713E+81 
3.5263E+8! 
3,9789E+8 | 
3.6292E+6 1 
3.6772E+81 
Seycoletel 
3.7669E+81 
3,8888E+83 
3.8498E+8! 


me eee 


2, 4957E+09 
2.9112E+88 
3.1478E+98 
3.3937E+88 
3.65) 26+83 
3.9219E+88 
4, 20435408 
4,4993E+88 
4, 8869E+88 
5.1266E+88 
5. 4588E+88 
5. 8802E+88 
6.1527E+88 
§.5146E+88 
6, 88546488 
7.2632E+88 
760306488 
8. 4495E+88 
8.4538E+88 
8.8657E+88 
9, 2848E+88 
9.7185E+88 
1,01426+81 
1.O579E+8! 
1. 1O21E+83 
1,1468E+81 
{,1919E+81 
1,2375€+61 
1, 28366481 
1,3302E+81 
1,3774E+81 
1,.4258E+81 
1.4736E+81 
1.5215E+81 
{.5783E+61 
1.4193E+81 
1,6685E+61 
1.7177E+81 
1,7669E+81 
1. 8159E+61 
{ .8647E+41 
Per isle+ ail 
1.961 1E+61 
2. 8885E+81 
2.8599E+81 
2. 1819E+81 
2.1477E+81 


173 
Ma 
177 
179 
181 
133 
185 
187 
189 
19] 
193 
199 
7 
199 
281 
283 
289 
287 
289 
211 
213 
219 
“7 
219 
221 
223 
229 
227 
229 
231 
233 
Zod 
237 
239 
241 
243 
249 
247 
249 
291 
pia 
290 
237 
299 
261 
263 
269 
267 
269 
271 
273 


8. 6888E-84 
2, 7B88E-84 
8, 88885 -84 
3. 9886E-84 
9, B880E-34 
9,1888E-84 
9, 2888E-84 
9, 3888E-84 


-1,{468E+88 
-1,1572E+88 
-{.1679E+88 
-{.17856+88 
-1.1892E+88 
-1.1993E+88 
~1.2898E+88 
~1,2193E+88 


9. 4888E-84 
9 O888E-84 


“1. 22915+88 
-1,.2383E+88 
9.6808E-34 -1,24736+88 
9, 788BE-84 -1.2564E+88 
9, 8868E-84 -1,2655E+98 
9. 9BBBE-84 -1.2737E+88 
1, 8808E-83 -1.2818E+88 
1.8188E-83 -1.2898E+88 
1.d288E-83 -1,.2982E+38 
1 .8380£-83 -1.3857E+88 
1. 8488E-83 -1.3129E+86 
1.8588E-83 -1.3208E+88 
1.8688E-83 -1.3272E+88 
1,8788£-83 -1.3338E+88 
1 .G8G0E-83 -1.3379E+88 
1. 8988E-83 -1.3468E+88 
1, 1888E-83 -1.35206+88 
1.1188E-83 -1,3579E+88 
1,1280£-83 -1.3633E+80 
1.1388E-83 -1,.3683E+83 
1.1488E-83 -1.3738E+88 
1.1988E-83 -1.3777E+88 
1.1680E-83 -1.3828E+88 
1.1788£-83 -1.3868E+88 
1,1888E-83 -1.3788E+88 
1.1988E-83 -1,3996E+88 
1.2828E-83 -1.3969E+88 
1. 2188E-83 -1.4008E+88 
1.22886-83 -1.4038E+88 
1.2380E-83 -1.4858E+88 
1,2408E-83 -1.4888E+88 
1,2508E-83 -1.4188E+88 
1. 2688E-83 -1.4119E+88 
1.2788E-83 -1.4137E+88 
1.2888E-83 -1.4151E+88 
1,.2988E-83 -1.4168E+88 
1, 3888E-83 -1.4166+88 
1.3188E-83 -1.4174E+88 
1.3288E-83 -1.4178E+88 
1,3388E-83 -1.4179E+88 
1,3488E-83 -1.41756+88 
1, 3588E-83 -1.41468E+88 
1,3688E-83 -1.41463E+88 


3.679048} 
3.7480E+81 
3. 8888E+8! 
3.8613E+8! 
Seo) se a 
3,9821E+81 
4. 8438E+6! 
4,1832E+81 
4, 14641648) 
4 ,2248E+81 
4, 2849£+81 
4 ,3457E+8i 
4,4862E+81 
4,4661E+81 
4,5263E+385 
4 ,58786+81 
4,46477E+81 
4. 7079E+61 
4.7676E+81 
4.8279£+81 
4, 8885E+81 
4, 9490E+81 
2. BO89E+81 
3.8684E+8 1 
v1 285E+81 
3.1889E+81 
3.2478E+81 
9. 3885E+8 1 
3.3676E+B1 
9,4274E+81 
5.4876E+81 
9.9473E+81 
3.6866E+81 
9. 6694E+81 
3. /248E+81 
9.7844E +81 
9.8437E+81 
9. 9825E+81 
3. 9606E+81 
6,8192E+81 
6,87815+8! 
6.1368E+81 
6.1958E+81 
6. 2524E +81 
6.3181E+8 
6,3686E+81 
6.42686 +81 
6.4846E+81 
6,9419E+81 
6.5998E+81 
6.46968E+8 1 


184 


1.360648} 
1,.3294E+81 
1. 2985E+8! 
1.24382E+61 
1.23885 +3; 
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